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PREFACE 


This  volume  describes  analyses  that  were  conducted  in  connection  with 
the  investigation  of  uncertainty  and  sensitivity  of  physiologically  based 
pharmacokinetic  (PP^K)  models.  In  particular,  Interest  was  focussed  on  the 
impact  of  uncertainty  about  the  values  of  u.odel  parameters  and  the  sensitivity 
of  PBPK  models  and  dose  response  models  to  the  values  of  the  parameters. 

Part  1  of  this  volume  describes  a  sensitivity  analysis  conducted  using  a 
relatively  simple  PBPK  model.  Parameter  values  were  varied  one  at  a  time  by 
an  arbitrary,  small  percentage.  The  percentage  change  in  PBPK  model  output 
(dose  surrogate  estimates)  was  recorded  for  each  parameter  change. 

Also  presented  in  Part  1  are  the  preliminary  considerations  and 
experimental  data  relevant  to  the  conduct  of  an  uncertainty  analysis.  In  such 
an  analysis,  the  parameters  are  allowed  to  vary  in  biologically  and 
experimentally  meaningful  ways  to  a  degree  consistent  with  the  observed 
uncertainty  and  variability  associated  with  the  parameter  values.  Then,  a 
distribution  of  output  values  (dose  surrogate  estimates  or  risk  estimates) 
derived  from  those  varying  parameters  is  available  for  subsequent  use.  The 
output  distribution  reflects  all  the  parameter  uncertainties. 

Part  2  presents  some  additional  analyses  related  to  the  uncertainty 
analysis  just  described.  The  contribution  of  an  Individual  parameter  (or  set 
of  interrelated  parameters)  to  the  overall  output  uncertainty  is  examined  by 
allowing  only  that  parameter  (or  set  of  parameters)  to  vary.  The  variation  In 
the  output  values,  relative  to  the  variation  in  the  output  when  all  parameters 
are  allowed  to  vary,  is  taken  as  an  indicator  of  the  contribution  of  that 
parameter  (or  set  of  parameters). 
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VOLUME  VI 


SENSmVITYAJNCERTAINTY  ANALYSIS 

PART  1  OF  2  PARTS 

SENSITIVITY  ANALYSIS; 


PRELIMINARY  CONSIDERATIONS  FOR  UNCERTAINTY  ANALYSES 


A.  INTRODUCTION 


An  analysis  of  PBPK  models  and  their  use  in  risk  assessment  should 
attempt  to  characterize  the  uncertainties  that  are  associated  with  that  use. 
Some  uncerta int ’ es  described  in  the  Introduction  to  this  document  (Volume  I, 
Part  1)  relate  to  the  more  or  less  general,  conceptual  problems  that  must  be 
resolved  in  order  to  make  the  use  of  PBPK  modeling  in  risk  assessment  less 
uncertain.  This  vol...  e  describes  an  investigation  of  a  more  specific  source 
of  uncertainty:  the  uncertainty  associated  with  the  estimation  of  parameters 
for  the  PBPK  model  and  the  sensitivity  of  the  desired  dose  surrogates  to 
changes  in  t  .ose  parameters.  The  investigation  focussed  on  the  simple  PERC 
model  discussed  in  Volume  III,  x ai t  1.  The  goal  was  to  relate  the  extent  of 
knowledge  about  the  input  parameters  to  the  range  of  dose  surrogate  values, 
and  ulti...ately  risk  estimates,  that  are  consistent  with  that  Knowledge. 

The  investigation  of  PBPK  model  uncertainty  was  divided  into  tvo  phases, 
a  sensitivity  analysis  and  an  uncertainty  analysis.  The  former  was  concerned 
with  the  degree  to  which  the  me iel  results  (dose  surrogate  estimates)  depended 
on  the  values  of  the  input  parameters  and  the  extent  of  change  in  those 
estimates  associated  with  changes  in  the  oarameter  values.  The  uncertainty 
analysis  examined  uncertainty  in  the  parameter  values  per  se  and  the 
associated  distribution  of  dose  surrogate  values  that  resulted  from  the 
simultaneous  consideration  of  all  the  uncertainties.  Of  course,  the  resulting 
distribution  of  the  dose  surrogate  values  depended  not  only  on  the  degree  of 
uncertainty  associated  with  a  parameter  but  also  on  the  sensitivity  of  the 
model  to  changes  in  that  parameter.  Thus,  the  sensitivity  analysis  was  the 
basis  for  '  oth  phases  of  this  investigation  and  is  described  first. 
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B.  SENSITIVITY 


The  sensitivity  analysis  consisted  of  an  examination  of  the  sensitivity 
of  the  PERC  model  to  changes  In  single  parameter  values.  Three  dose 
surrogates  were  used  to  measure  sensitivity:  area  under  the  liver 
concentration  curve  of  the  parent  (AUCL) ,  area  under  the  arterial  blood 
concentration  curve  of  the  parent  (AUCA) ,  and  virtual  concentration  of  the 
metabolite  in  the  liver.  The  average  daily  values  of  these  dose  surrogates 
were  used.  Three  sets  of  "preferred"  parameter  values  were  selected,  one  for 
mice,  one  for  rats,  and  one  for  humans  (Table  VI-1-1).  In  addition,  because 
it  is  likely  that  route  of  exposure,  pattern,  and  level  of  exposure  may  modify 
the  sensitivity  of  the  model  to  a  particular  parameter,  dosing  variations  were 
defined  by  combining  several  options  with  respect  to  such  aspects.  The  routes 
and  patterns  used  were; 

•  inhalation,  exposure  for  8  hours  per  day 

•  inhalation,  continuous  exposure 

•  gavage ,  once  per  day 

•  intravenous  Injection,  once  per  day. 

For  each  of  the  four  route/pattem  combinations  three  separate  dose  levels 
\/ere  chosen.  A  low  level  of  exposure  was  selected  such  that  the  metabolic 
pathway  was  far  from  saturation.  An  intermediate  dose  level  at  which  the 
metabolic  pathway  was  changing  from  first-order  to  zero-order  and  a  saturating 
dose  provided  the  other  two  doses.  Preliminary  runs  of  the  model  with  the 
parameter  sets  listed  in  Table  VI -1-1  allowed  estimation  of  the  appropriate 
dose  levels  (Table  VI-1-2). 
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For  the  gavage ,  Intravenous,  and  8 -hour  inhalation  exposure 
scenarios,  the  model  was  run  to  simulate  24  hours  starting  at  the  beginning  of 


exposure.  For  the  continuous  inhalation  scenario,  an  initial  time  was 
determined  such  that  at  that  time  the  system  was  close  to  steady- state .  The 
model  was  then  allowed  to  simulate  another  24  hours.  It  was  that  24-hour 
period  over  which  the  PERC  liver  concentration  and  the  arterial  blood 
concentration  were  integrated  and  over  which  the  metabolism  was  summed  to 
yield  the  dose  surrogates  of  interest.  The  times  at  which  the  appropriate 
periods  began  are  listed  in  Table  VI-1-2. 

For  each  of  the  sets  of  "preferred"  parameter  values  and  for  each 
selected  combination  of  route,  pattern,  and  level  of  dosing,  the  model  was  run 
to  obtain  baseline  values  for  each  of  the  dose  surrogates.  Then,  each 
parameter  (one  at  a  time)  was  increased  by  1%:  when  examining  the  sensitivity 
to  parameter  X,  all  other  parameters  remained  fixed  at  their  "preferred" 
values.  A  minor  exception  to  this  rule  related  to  the  handling  of  blood  flow 
rates.  The  flows  to  the  tissue  compartments  must  sum  to  the  total  cardiac 
output.  Thus  when  one  compartment  flow  rate  was  increased  by  1%,  the  total 
caraiac  output  was  increased  by  the  same  amount.  The  other  individual 
compartment  flow  rates  were  not  changed.  Thus  the  effect  of  this  treatment 
was  to  change  slightly  the  percentage  of  the  cardiac  output  that  reaches  any 
particular  compartment. 

Sensitivity  was  expressed  in  terms  of  the  percent  change  in  the  values 
of  the  three  dose  surrogates.  If  DS  was  one  of  the  dose  surrogates,  DSg  was 
its  baseline  value,  and  DSj  wa.s  Its  value  obtained  by  increasing  the  value  of 
parameter  X  by  1%,  then  100*(DS,  -  DSo)/DSo  was  the  percent  change  recorded. 
Tables  VI  - 1  -  3  through  V-1-14  present  the  results  of  the  sensitivity 
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investigation  of  the  PERC  model.  [In  those  tables,  percent  changes  less  than 
10'^  have  been  uniformly  designated  by  asterisks  to  indicate  changes  at  or 
below  the  accuracy  of  the  model.] 

Some  patterns  were  evident  from  examination  of  the  values  reported  in 
the  tables.  Not  surprisingly,  the  AUCL  dose  surrogate  was  most  sensitive  to 
pi  for  almost  all  routes  of  exposure,  dose  levels,  or  species.  The 
liver/blood  partition  coefficient,  pi,  determines  the  extent  to  which  PERC 
concentrates  in  the  liver.  For  some  inhalation  exposures  at  high  enough  dose 
(Table  Vl-1-9),  the  sensitivity  of  AUCL  to  pb  exceeded  that  of  AUCL  to  pi. 

This  must  be  because  the  blood  becomes  so  laden  with  PERC  that  a  large  amount 
of  the  chemical  will  partition  to  the  liver  for  both  pi  values  examined.  For 
intravenous  administration,  AUCL  was  more  sensitive  to  the  parameter  tlv  (the 
duration  of  an  iv  dose  administration)  than  to  pi  in  rats  and  humans  at  medium 
doses  (Tables  VI-1-8  and  VI-1-12). 

For  AUCA,  pb  and  qpc  (which,  in  general,  were  very  similar  for  gavage  or 
intravenous  dosing  in  the  magnitude  of  the  sensitivity  values  but  differed  in 
sign)  were  the  most  important  parameters.  [Note:  in  this  discussion, 
parameters  to  which  a  dose  surrogate  is  more  highly  sensitive  are  referred  to 
as  "important"  parameters.  To  a  considerable  degree  this  designation  is 
relative,  depending  as  it  does  on  the  sensitivity  of  the  dose  surrogate  to 
other  parameters.)  For  virtual  concentration  of  metabolite  (CM),  vie  and, 
especially  for  the  human  parameter  set,  vmaxc  were  most  important. 

The  body  weight  parameter  (bw)  displayed  an  Interesting  pattern.  AUCL 
and  AUCA  were  more  sensitive  to  bw  than  was  CM  when  dose  was  given  by  gavage 
or  iv  but  these  two  dose  surrogates  were  less  sensitive  than  was  CM  when 
inhalation  exposures  were  considered.  This  pattern  was  especially  evident  for 
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the  mouse  parameter  set.  It  is  suspected  that  this  pattern  was  due,  at  least 
in  part,  to  the  manner  in  which  doses  were  expressed.  The  Iv  and  gavage 
scenarios  were  defined  in  terms  of  doses  given  in  mg/kg  body  weight  Thus, 
changes  in  body  weight  also  changed  the  dose  adjuinistered  by  these  two  routes 
but  not  for  the  inhalation  exposures.  The  sensitivity  of  CM  in  the  latter 
cases  is  most  likely  a  reflection  of  the  differential  scaling  of  vmaxc  (scaled 
according  to  bw^'  and  vl  (scaled  according  to  bw^  °)  .  Since  all  the  flows 
and  compartment  volumes  are  scaled  proportionally  to  bw^°,  the  effects  on 
blood  and  liver  PERC  concentrations  due  to  a  change  in  body  weight  were  not  as 
substantial . 

As  mentioned  above,  vmaxc  was  important  for  CM.  It  was  substantially 
less  important  to  AUCL  and  AUCA,  although  for  mice,  particularly  at  low  doses, 
its  effect  on  these  two  dose  surrogates  can  be  relatively  important  and  even 
exceeded  the  effect  on  CM  (Tables  VI-1-3  through  VI-1-6).  It  was  also  the 
case  that  the  sensitivity  of  CM  to  vmaxc  increased  as  dose  Increased. 

However,  the  sensitivity  of  AUCA  and  AUCL  to  vmaxc  decreased  as  dose  level 
increased.  As  doses  increased  and  as  metabolism  became  saturated,  entailing 
relatively  less  PERC  being  metabolized,  the  effect  of  metabolism  on  the  parent 
concentrations  became  less  important.  Note  that,  of  the  three  parameter  sets 
considered,  the  mouse  set  had  the  largest  vmaxc,  followed  by  the  rat  set.  As 
vmaxc  decreased,  the  sensitivity  of  CM  to  vmaxc  increased  while  the 
sensitivity  of  AUCL  and  AUCA  to  vmaxc  decreased. 

The  other  metabolic  parameter,  km,  became  less  important  as  dose  level 
increased,  as  expected.  This  was  especially  true  for  the  dose  surrogates  AUCL 
and  AUCA.  Those  two  dose  surrogates  were  more  sensitive  to  km  for  mice  than 
for  rats  or  humans.  On  the  other  hand,  at  low  to  medium  dose  levels,  CM  was 
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less  sensitive  to  km  when  the  parameters  were  equal  to  the  values  in  the  mouse 
set  than  when  they  took  the  rat  or  human  values.  The  km  value  for  rats  and 
humans  was  about  two  times  larger  than  the  value  for  mice. 

All  the  dose  surrogates  examined  were  relatively  insensitive  to  changes 
in  the  partition  coefficients  corresponding  to  compartments  other  than  the 
liver,  regardless  of  route  or  level  of  exposure.  In  fact,  pi  was  Important 
only  for  the  AUCL  surrogate.  Conversely,  pb  was  relatively  important  for  all 
routes,  species,  and  dose  levels.  The  sensitivity  of  the  dose  surrogates  to 
pb  tended  to  increase  with  dose  level,  although  the  differences  across  dose 
levels  were  generally  not  great.  An  exception  to  this  observation  was  noted 
for  the  CM  surrogate,  for  which  the  sensitivity  to  pb  tended  to  decrease  at 
high  doses,  especially  for  continuous  inhalation  exposures.  In  those  cases, 
the  high  doses  entailed  high  enough  blood  concentrations  that  the  limiting 
step  was  metabolism,  not  delivery  of  PERC  to  the  liver. 

For  a  compartment  other  than  the  liver,  the  sensitivity  of  the  dose 
surrogates  to  the  partition  coefficient  and  to  the  compartment  volume  were 
almost  indistinguishable.  This  relationship  was  strongest  at  the  high  doses 
but  was  evident  also  at  medium  and,  to  a  lesser  extent,  low  dose  levels. 

(This  observation  may  have  important  implications  for  the  optimization  of 
parameters  in  the  face  of  uncertainty  about  both  compartment  volumes  and 
partition  coefficients.  Generally,  it  will  be  the  case  that  a  partition 
coefficient  value  is  more  uncertain  than  the  corresponding  compartment 
volume.)  The  observation  that  partition  coefficients  other  than  pb  and  pi 
were  relatively  unimportant  extends  to  compartment  volumes  as  well.  Indeed, 
dose  surrogates  other  than  CM  were  not  very  sensitive  to  liver  volume  changes. 
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certainly  no  more  so  than  to  changes  ^n  the  other  compartment  volumes.  The 
importance  of  vie  to  the  CM  surrogate  was  noted  above. 

The  ventilation  rate  parameter  was  relatively  Important  for  all  dose 
surrogates  and  all  routes  of  exposure  at  low  to  medium  dose  levels.  The 
sensitivity  of  the  dose  surrogates  to  ventilation  rate  (qpc)  varied  somewhat 
with  the  level  of  dose.  Moreover,  the  dose  trend  that  was  observed  depended 
on  the  route  of  exposure.  For  Inhalation  exposures,  the  surrogates  became 
less  sensitive  to  qpc  as  dose  increased  (at  high  enough  doses  the 
concentration  in  the  blood  was  high  no  matter  how  fast  the  animals  took  in  the 
chemical).  For  gavage  and  intravenous  dosing,  the  sensitivity  increased  as 
dose  level  increased  (high  iv  or  gavage  doses  saturated  the  bIcod  and 
exhalation  was  a  prime  means  of  eliminating  the  parent  but  not,  as  in  the  case 
of  inhalation,  of  taking  in  PERC) .  As  might  be  expected,  increasing  qpc 
increased  the  dose  surrogate  estimates  when  dose  was  administered  via 
inhalation  whereas  an  increase  in  qpc  reduced  the  surrogate  estimates  when 
breathing  rate  did  not  influence  the  absorption  of  the  chemical,  i.e.,  when 
PERC  was  given  via  gavage  or  iv. 

As  noted  above,  qpc  was  a  more  important  parameter  for  AUCA  than  for  the 
other  surrogates.  This  was  not  due  to  qpc  sensitivity  values  that  were 
substantially  greater  for  AUCA  than  for  the  other  surrogates.  Rather,  the 
other  parameters  that  were  important  for  AUCL  and  CM  (pi,  vie,  and  vmaxc,  for 
example)  were  much  less  important  for  AUCA  so  that  the  relative  Importance  of 
qpc  to  AUCA  increased. 

The  blood  flow  rates  were  among  the  least  Important  parameters  for  any 
dose  surrogate.  It  is  interesting  to  note  the  changes  in  the  direction  of 
change  of  a  dose  surrogate  Induced  by  an  Increase  in  a  parameter  value.  For 
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Increases  in  the  blood  flow  rates,  it  was  frequently  the  case  that  the  signs 
on  the  sensitivity  values  differed  for  different  dose  levels,  parameter  sets, 
and  dose  surrogates.  No  consistent  pattern  was  apparent.  This  effect  was  not 
anticipated. 

Another  way  to  axamlne  the  sensitivity  issue  is  to  rank-order  the 
parameters  with  respect  to  the  degree  of  sensitivity  of  a  dose  surrogate  to 
those  parameters.  Tables  VI-1-15  and  VI-1-16  display  the  orderings  obtained 
for  the  CM  surrogate  when  two  routes  of  exposure  and  two  species  (parameter 
sets)  were  considered.  The  importance  of  vie  and  vmaxc  to  CM  is  underscored 
in  these  representations.  Several  other  observations  are  relevant. 

•  The  Interaction  of  dose  route  and  dose  level  in  defining 
sensitivities  was  seen,  for  example,  in  the  increasing  importance  of 
qpc  as  iv  dose  increased  and  as  inhalation  dose  decreased  for  the  rat 
parameter  set. 

•  The  distinction  between  important  and  unimportant  parameters  was  much 
clearer  whe  \  a  high  inhalation  dose  was  administered,  whereas  a  more 
continuous  variation  in  sensitivity  was  seen  at  lower  doses  or  for 
other  routes. 

•  Blood  flow  rates  to  the  compartments  were  consistently  unimportant 
(with  the  possible  exception  of  qlc)  . 

There  was  little  difference  between  the  two  parameter  sets  with  respect 
to  the  most  sensitive  parameters,  although  the  particular  ordering  differed 
somewhat  between  the  species.  This  was  a  reflection  of  the  difference  in  the 
values  assigned  to  the  parameters.  It  was  not  the  case,  however,  that  the 
sole  determinant  of  the  importance  of  a  parameter  was  the  value  of  that 
parameter.  Consider  the  parameters  defining  the  percentage  of  cardiac  output 
directed  to  the  compartments.  Even  though  these  parameters  were  the  same  in 
the  rat  and  mouse  parameter  sets,  the  ranks  that  these  parameters  received 
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differed  across  species.  So  too  did  the  ranking  based  solely  on  those  blood 
flow  parameters,  with  qsc  ranked  higher  than  qfc  In  one  species  or  for  one 
dose  level  and  vice  versa  for  the  other  parameter  set  and  for  other  dose 
levels . 

As  stated  above,  the  sensitivity  of  the  dose  surrogates  to  the 
parameters  is  an  important  factor  in  the  determination  of  uncertainty 
propagation.  Those  parameters  for  which  uncertainty  with  respect  to  their 
real  value  is  great  may  not  necessarily  lead  to  great  uncertainty  in  the  dose 
surrogates  if  the  surrogates  are  not  particularly  sensitive  to  those 
parameters.  On  the  other  hand,  even  relatively  minor  uncertainty  regarding 
the  value  of  a  parameter  may  entail  large  uncertainty  in  a  dose  surrogate  (and 
ultimately  in  risk  estimates  based  on  dose  surrogate  estimates)  if  the 
surrogates  of  interest  are  extremely  sensitive  to  the  parameter  in  question. 
Uncertainty  is  discussed  more  completely  in  the  next  section. 

C .  UNCERTAINTY 


The  values  of  the  parameters  that  define  a  PBPK  model  are  uncertain. 

The  uncertainty  associated  with  a  parameter  estimate  may  be  due  to  lack  of 
adequate  data,  differences  in  methods  used  to  estimate  the  parameter  value, 
and/or  recognition  of  the  fact  that  a  single  value  may  not  be  adequate  to 
characterize  a  parameter  that  varies  over  the  population  of  interest.  The 
goal  of  the  uncertainty  analysis  described  here  was  to  relate  the  parameter 
uncertainties  to  the  distribution  of  the  results  of  the  PBPK  models,  l.e.,  to 
the  dose  surrogates  estimated  by  those  models.  In  addition,  that  analysis 
Included  computation  of  distributions  of  risk  estimates,  which  depended  on  the 
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values  estimated  for  the  dose  surrogates,  and  thus  were  also  related  to 
parameter  uncertainty.  We  refer  to  the  distribution  of  dose  surrogates 
attributable  to  parameter  uncertainty  as  dose  surrogate  uncertainty. 

Similarly,  the  distribution  of  risk  estimates  attributable  to  parameter 
uncertainty  is  referred  to  as  risk  uncertainty.  The  definition  of  risk 
uncertainty  is  extended  to  Include  effects  of  uncertainty  about  the  true  rate 
of  response  for  any  group  of  animals  tested  (observed  for  tumor  response)  and 
used  for  fitting  of  a  dose  -  response  model. 

The  first  task  in  this  uncertainty  analysis  was  the  definition  of  the 
joint  uncertainty  distribution  for  the  parameters.  The  description  of  that 
work  for  the  PERC  model  parameters  comprises  the  bulk  of  this  section.  As 
described  below,  the  parameter  uncertainty  distribution  takes  into  account  as 
many  of  the  features  of  the  estimation  of  parameter  values  as  possible.  The 
emphasis  in  this  section  is  on  development  of  methods  and  conceptual 
approaches  to  the  issue  of  parameter  uncertainty  estimation. 

Given  estimated  uncertainty  distributions  for  the  parameters,  the  dose 
surrogate  uncertainty  distributions  were  estimated  by  Monte  Carlo  simulation. 
For  each  simulation,  a  value  was  sampled  from  each  of  the  distributions 
defined  for  the  parameters.  These  values  were  then  used  in  the  PBPK  model  to 
determine  the  corresponding  values  of  the  dose  surrogates.  A  large  number  of 
simulations  allowed  estimation  of  the  uncertainty  distribution  of  the  dose 
surrogates . 

Dose  surrogate  uncertainty  was  then  extended  to  estimate  the  uncertainty 
in  risk  estimates  due  to  uncertainties  in  model  parameters.  The  uncertainty 
distributions  for  the  animal  surrogate  doses  were  estimated  as  described  in 
the  preceding  paragraph,  for  each  of  the  dose  patterns  used  in  one  or  more 
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carcinogenicity  bioassay.  Similarly,  uncertainty  distributions  for  human 
surrogate  doses  were  estimated  for  (the  usually  low)  levels  of  exposure  that 
are  of  concern;  this  was  accomplished  by  sampling  from  the  parameter 
uncertainty  distributions  corresponding  to  the  human  parameter  values  and 
running  the  PBPK  model  many  times  exactly  as  described  above  for  the  animal 
parameters  and  simulations. 

Given  the  aistributions  of  the  dose  surrogate  estimates  (for  each  animal 
dosing  pattern  and  for  the  human  exposure  scenario  of  Interest),  it  was 
possible  to  estimate  risk  uncertainty.  The  methods  by  which  this  was 
accomplished  for  one  carcinogenicity  bioassay  are  described  as  follows. 


1.  A  bioassay  was  selected.  Carcinogenic  responses  appropriate  for 
use  in  estimation  of  human  risk  were  chosen. 

2.  For  the  doses  employed  in  the  bioassay,  the  dose  surrogate  values 
co’'responding  to  the  preferred  PBPK  model  parameter  values  were 
used  to  fit  a  dose  -  response  curve. 

3.  In  order  to  incorporate  the  uncertainty  about  the  true  response 
rate  at  the  doses  tested,  the  probabilities  of  response  estimated 
by  the  dose  -  response  curve  fit  In  step  2  were  used  to  randomly 
regenerate  response  rates  (from  binomial  distributions  with  their 
parameters  defined  by  the  number  of  animals  tested  at  each  dose 
and  the  estimated  probabilities  of  response). 

4.  The  dose  surrogate  distributions  (previously  estimated)  were 
sampled  once  for  each  dose  group  and  a  new  dose-response  function 
was  fit  using  the  just-sampled  dose  surrogate  values  and  the 
randomly  generated  tumor  response  rates. 

5.  A  value  from  the  uncertainty  distribution  of  the  human  dose 
surrogates  (also  determined  oreviously)  was  sampled;  it  determined 
the  value  for  which  risk  was  estimated,  which  was  accomplished 
using  the  newly  fitted  dose-response  curve. 


Repetition  of  this  procedure  a  large  number  of  times  allowed  estimation  of  the 
uncertainty  distribution  of  the  human  risk  estimates. 
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Several  results  are  available  from  an  uncertainty  analysis  as  described: 

•  Uncertainty  distributions  for  each  of  the  parameters:  the  uncertainty 
distribution  for  each  parameter  or  set  of  parameters  is  represented  by 
a  statistical  probability  distribution.  Several  of  the  parameters  for 
the  PERC  model  are  represented  in  this  manner  below. 

•  Uncertainty  distributions  for  dose  surrogates:  the  empirical 
distributions  resulting  from  the  Monte  Carlo  sampling  of  the  parameter 
uncertainty  distributions  can  be  displayed  graphically  as  cumulative 
frequency  plots,  histograms,  or  displayed  in  tabular  form,  listing, 
for  example,  the  2.5,  5,  10,  90,  95,  and  97.5  percentiles. 

•  Uncertainty  distributions  for  the  risk  estimates:  the  empirical 
distributions  resulting  from  Monte  Carlo  sampling  of  the  parameter 
uncertainty  distribution,  PBPK  modeling  for  each  bioassay  dose  level 
and  human  exposure  level,  and  dose  -  response  analysis. 


The  remainder  of  this  section  reviews  the  considerations  and  data  relevant  to 
estimation  of  uncertainty  distributions  for  the  parameters  of  the  PERC  PbPK 
model  considered  in  Volume  III,  Part  1  (and  used  for  the  sensitivity  analysis 
in  this  volume) .  The  focus  has  been  on  the  development  of  the  methodological 
approaches  and  ider.'-  fication  of  difficulties  that  may  arise  in  the 
specification  of  uncertainty  distributions. 

Appendix  VI-l-A  contains  a  published  analysis  that  has  completed  the 
task  of  deriving  dose  surrogate  and  risk  estimate  uncertainty  distributions. 
The  example  in  that  document  is  also  for  PERC,  and  uses  results  from  a  female 
mouse  carcinogenicity  bloassay  to  get  risk  estimates. 


1_, _ Parameter  Uncertainty  Estimation 

The  uncertainty  associated  with  a  parameter  value  depends  on  the  type  of 
data  available  for  estimating  the  parameter  and  the  way  that  the  data  are  used 
to  arrive  at  an  estimate.  Frequently  more  complex  and  indirect  estimates 
(e.g.,  estimates  based  on  numerous  assumptions)  will  be  attended  by  greater 
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uncertainty  than  relatively  simple  and  direct  estimates.  For  example,  if  one 
attributes  to  humans  a  value  measured  directly  only  in  rats,  then  the 
uncertainty  for  the  hviman  value  is  affected  by  the  degree  of  variation  in  the 
value  from  species  to  species,  as  well  as  whatever  uncertainties  are 
associated  with  the  direct  measurement  in  rats. 

In  an  uncertainty  analysis  as  described  here,  it  is  necessary  to  have 
"preferred"  values  for  parameters.  In  general,  a  preferred  value  for  a 
parameter  is  an  estimate  that  has  been  determined  by  one  or  more  users  of  the 
specific  PBPK  model  and  has  been  validated  to  some  extent.  The  preferred 
value  roughly  corresponds  to  the  center  of  the  uncertainty  distribution  that 
is  to  be  estimated.  For  the  PERC  model  parameters,  for  example,  the  preferred 
values  may  be  assumed  to  be  those  given  by  Reitz  and  Nolan  (1986)  (cf.  Table 
III-l-l)  . 

The  manner  in  which  Reitz  and  Nolan  derived  estimates  for  the  parameters 
of  their  PERC  model  is  representative  of  the  considerations  that  typically 
apply  to  the  issue  of  parameter  estimation.  The  necessity  of  attributing  to 
humans  parameter  values  based  on  measurements  in  animals  is  but  one  example. 
Thus,  the  approaches  that  they  have  used  are  convenient  starting  points  for 
this  investigation  of  parameter  uncertainty  estimation  and  are  examined  in 
some  detail  here.  The  approaches  of  Reitz  and  Nolan  do  not  necessarily 
provide  the  most  reliable  estimates  possible  for  PBPK  model  parameters; 
however,  by  examining  their  approaches,  one  can  evaluate  the  uncertainties 
that  will  arise  in  many  situations.  One  nice  feature  of  the  approach  to 
uncertainty  estimation  described  below  is  that,  as  more  reliable  estimates  of 
parameter  values  become  available,  the  uncertainties  associated  with  the 
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parameteii;  will  be  reauced  because  of  the  ccpslderations  given  to  the  sources 
for  and  assiunptions  inherent  in  the  estimation  procedure. 

Probaoilltv  Models  for  Uncertainty.  The  first  assumption  underlying 
this  approach  to  uncertainty  analysis  is  that  uncertainties  are  naturally 
expressed  in  terms  or  probability.  For  example,  if  one  is  fairly  certain  that 
a  parameter  value  must  lie  within  a  specified  inteir'al  then  that  interval  will 
be  assumed  to  contain  the  true  parameter  value  with  high  probability.  (Note; 
for  continuous  variables  the  probability  of  any  single  value  is  zero,  while 
the  probabl’ Ity  density  of  the  value  may  be  greater  than  zero.  A  probability 
density  may  be  integrated  over  an  interval  to  yield  the  probability  associated 
with  the  interval.)  ("iven  this  basic  assumption,  the  formal  treatment  of 
uncertainties  involves  specification  of  a  joint  probability  distribution  for 
all  ^he  parameters  tbit  are  subject  to  uncertainty. 

Consider  an  example  for  a  p.-'rameter  that  is  assumed  to  be  independent  of 
the  other  parameters.  A  parameter  is  independent  of  another  para’'  ,r  when 
Information  regarding  the  value  of  the  other  parameter  is  uninformative  with 
respect  to  the  value  of  the  first  parameter.  In  this  case,  the  preferred 
value  of  the  parameter  might  be  assigned  the  highest  probability  density,  and 
values  at  Increasing  distances  (both  above  and  below  the  preferred  value)  are 
associated  with  lower  density.  If  it  is  further  assumed  that  valu  s 
equidistant  from  the  preferred  value  have  equal  probability  density  and  that 
every  value  has  nonzero  (positive)  density,  then  a  natural  representation  of 
the  uncertainty  distribution  is  given  by  the  normal  distribution.  The  mean  of 
the  distribution  is  equated  with  the  preferred  value  and  the  variance  (the 
only  other  parameter  defining  a  normal  curve)  is  estimated  on  the  ba.sis  of  tne 
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data  used  to  deii’'’e  the  parameter  estimates.  The  variance  would  be  the 
probability  parameter  reflecting  uncertainty. 

The  probability  distribution  for  a  single  parameter  is  called  the 
marginal  probability  distribution  ror  that  parameter.  A  set  of  marginal 
probability  distributions  (one  for  each  parameter)  specifies  the  joint 
probability  distribution  when  the  parameters  are  independent.  Wheii  the 
parameters  are  i.ot  assumed  to  be  independent,  additional  assumptions  or 
parameters  may  be  required  to  specify  the  nature  of  the  dependencies  and  to 
characterize  the  joint  distribution. 

The  assumptions  discussed  above  in  relation  to  the  representation  of  the 
uncertainty  distribution  of  a  nararaeter  by  a  normal  curve  are  generally  not 
appropriate  for  characterization  of  PBPK  model  parameter  uncertainty.  This  is 
true  because  of  the  known,  '’ogical  bounds  on  the  values  of  the  parameters. 
Obviously,  all  the  parameters  are  constrained  to  be  nonnegative  Moreover, 
certain  of  the  parameters  (such  as  percentage  of  body  weight  that  is  occupied 
by  a  compartment  and  proportions  of  cardiac  output  that  flow  to  a  compartment) 
are  bounded  above  by  unity  (or  100%).  The  parameters  that  are  proportions  or 
percentages  are  also  constrained  by  the  values  of  the  other  parameters  that 
describe  proportions  associated  with  other  compartments.  Some  distributions 
that  may  be  considered  for  defining  parameter  uncertainty  are  described  below; 
they  all  are  expressed  in  terms  of  ixplicit  probability  models  that  have  a 
finite  mean  (generally  equated  with  the  preferred  value  of  the  corresponding 
parameter)  and  addi,.ional  parameters  specifying  the  distribution  of 
probability  about  the  m.e<'in. 

Two  cases  arc  recognized.  For  parameters  that  are  logically  bounded 
below,  one  may  use  the  lognormal  probability  distribution,  l.e.,  one  assumes 
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that  uncertainty  follows  a  normal  distribution  with  a  specified  mean  anJ 
variance  when  parameter  values  are  expressed  in  the  log  scale.  The  lognormal 
is  flexible  in  the  sense  that  different  selections  of  log-scale  mean  and 
variance  may  lead  to  distributions  with  widely  different  shapes,  although  all 
assign  zero  probability  density  to  values  zero  or  smaller.  If  the  variance  is 
small  relative  to  the  mean,  then  the  shape  of  the  density  function  is  similar 
to  that  of  the  normal  distribution;  relatively  larger  values  of  the  variance 
lead  to  J -shaped  distributions. 

For  lognormal  distributions,  one  can  assume  that  the  mean  of  the 
distribution,  in  the  log  scale,  is  the  log  of  the  preferred  value.  Then  for  a 
preferred  value  X  and  any  positive  number  K,  the  values  X/K  and  X*K  have  equal 
probability  density  (i.e.,  the  probability  density  function  is  symmetric  in 
the  log  scale  about  the  log  of  the  preferred  value.)  It  is  convenient  to 
represent  uncertainty  using  an  "uncertainty  factor"  (UF)  not  smaller  than 
unity  iuch  that  for  the  preferred  value  X,  the  Interval  X/UF  to  X*UF  Is 
considered  to  contain  the  true  value  of  the  parameter  with  probability  at 
least  95%.  More  about  uncertainty  factors  and  their  estimation  is  given 
below. 

In  the  second  case,  for  parameters  that  are  proportions  (e.g. 
compartment  volumes  as  percentage  of  total  body  volume),  a  joint  distribution 
must  be  specified  for  each  set  of  proportions,  since  the  observed  proportions 
cannot  be  varied  independently.  One  can  use  a  Dirichlet  distribution  for  each 
set  of  proportions  (Johnson  and  Kotz,  1972),  which  is  a  multivariate 
generalization  of  the  beta  distribution.  (The  beta  distribution  is  the  most 
commonly  encountered  continuous  distribution  for  a  single  proportion:  Iman  and 
Shortencarier  (1984)  have  previously  discussed  the  use  of  the  beta 
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distribution  in  the  contexts  of  uncertainty  analysis.)  For  each  proportion  or 
percentage,  the  Dlrichlet  distribution  permits  only  values  between  0  and  1. 
Subject  to  this  constraint,  the  Dlrichlet  distribution  is  extremely  flexible 
in  terras  of  the  diversity  of  shapes  which  may  result  from  selection  of  various 
values  for  the  parameters.  Given  a  preferred  value  for  each  proportion,  the 
variances  and  covariances  of  the  proportions  are  determined  by  a  parameter 
THETA.  A  brief  discussion  of  the  Dlrichlet  distribution  is  provided  here. 

For  purposes  of  discussion,  assume  that  the  parameters  are  proportions, 
i.e.,  that  their  sum  is  constrained  to  equal  1.  For  K  proportions,  it  is 
necessary  to  specify  K  parameters  for  the  Dlrichlet  distribution.  There  are 
two  common  ways  of  expressing  the  parameters  of  a  Dlrichlet  distribution.  One 
representation  is  given  in  terms  of  the  parameters  p(l) ,p(2) , . . . ,p(K-l) 

(where  p(i)  is  the  expected  value  of  the  ith  proportion)  and  the  parameter 
THETA  (which  can  be  no  smaller  than  zero)  .  Note  that  the  Kth  expectation  is 
known  if  the  first  K-1  are  specified.  A  second  common  representation  is  in 
terms  of  a(i)  -  p(i)*THETA,  for  i  -  1,...,K.  The  first  representation  is  more 
natural  to  PBPK  modeling  applications.  The  preferred  value  for  the  ith 
proportion  ie  set  equal  to  p(i):  the  parameter  THETA  must  be  estimated  from 
experimental  data. 

It  is  clear  that  two  proportions  (such  as  proportion  of  body  weight  that 
is  the  liver  and  proportion  that  is  fat)  must  be  negatively  correlated.  That 
is,  as  the  proportion  of  body  weight  that  is  liver  Increases,  the  proportion 
that  is  fat  must  tend  to  decrease.  In  fact,  the  covariance  of  two  proportions 
described  Jointly  by  a  Dlrichlet  distribution  is  given  by 
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(1)  COV(i.J)  -  -  p(l)*p(J)/(THETA  +  1). 


Thus  the  uncertainties  for  members  of  a  set  of  proportions  assumed  to  be 
Dirichlet-distributed  are  not  independent.  They  are  completely  specified  by 
the  preferred  values  and  the  parameter  THETA. 

As  an  example,  consider  volumes  of  tissue  compartments.  These  can  be 
represented  by  a  Dirichlet-distribution  when  they  are  expressed  as  proportions 
of  the  total  body  weight.  In  fact,  the  proportions  may  change  with  changes  in 
body  weight.  For  instance,  it  appears  reasonable  that  the  volume  of  the  fat 
compartment  will  be  correlated  with  body  weight.  A  relatively  simple 
representation  of  the  dependence  of  fat  volume  on  body  weight  is  given  by 


(2a) 

pF(BW) 

-  pF*BW/BWp 

for  BW  <  BWp; 

(2b) 

pF(BW) 

-  1  -  (BWp/BW)*(l  -  pF) 

for  BW  1  BWp, 

where  pF  is  the  preferred  value  for  proportion  that  is  fat  at  the  preferred 
body  weight,  BWp.  This  formula  implies  that  the  expected  value  of  fat  volume 
goes  down  quadratically  with  body  weight  for  body  weight  less  than  BWp,  and 
that  body  weight  in  excess  of  BWp  consists  of  fat.  The  preferred  values  of 
the  remaining  (non-fat)  proportions,  p(l),  can  be  adjusted  by  multiplying  each 
by  (1 -pF(BW) )/( 1 -p( i) ) ,  so  that  compartment  proportions  other  than  that  for 
the  fat  compartment  retain  the  same  relationship  among  themselves  that  they 
have  at  the  preferred  body  weight.  The  pF(BW)  and  other  adjusted  p(i)'s  are 
the  parameters  that,  along  with  THETA,  now  define  the  Dirichlet  distribution 
for  compartment  proportions  in  an  animal  of  weight  BW. 
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In  summary,  then,  PBPK  model  parameter  uncertainty  can  be  modeled  on  the 
basis  of  probability  distributions.  For  parameters  other  than  those  that  are 
expressed  in  terms  of  proportions,  log-normal  distributions  can  be  assumed. 

For  proportions,  the  Dirichlet  distribution  provides  a  convenient  definition 
of  the  relationship  among  the  proportions.  As  illustrated  in  the  case  of 
compartment  volumes  as  proportions  of  body  weight,  inter- relations^'ips  among 
parameters  of  a  PB?K  model  can  be  incorporated  so  as  to  represent  biological 
real ities . 

What  remains  is  to  specify  tlie  relevant  sources  of  uncertainty  and  the 
manner  in  which  they  car.  be  used  to  estimate  the  other  parameters  of  the 
probability  distributions  (such  as  the  log-variance  or  leg-standard  deviation 
in  rhe  case  of  log-normal  distributions) ,  the  ones  that  actually  define  the 
uncertainty  about  preferred  values. 

Sources  of  Uncertainty.  Two  general  types  of  uncertainty  appear  to  be 
relevant  to  the  definition  of  the  uncertainty  distributions  discussed  above. 
The  first,  which  will  be  denoted  as  source  I,  is  uncertainty  about  the  value 
of  a  parameter  for  a  random  individual  in  a  population  that  is  due  to 
individual  variability  about  a  population  meaii.  This  uncertainty  may  be 
reflected  in  the  variation  in  paidmeter  estimates  among  individual  subjects 
within  a  study  and  expressed  in  terms  of  the  standard  deviation  or  coefficient 
of  variation.  The  second  source  of  uncertainty  is  identified  with  uncertainty 
regarding  the  population  mean  value  of  a  parameter;  this  source  is  denoted  as 
source  M.  Data  from  different  studies  that  suggest  different  mean  values  are 
relevant  to  this  type  of  uncertainty.  Both  sources  are  pertinent  to  the 
estimation  of  total  uncertainty  and  can  be  incorporated  by  considering  the 
following  framework. 
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The  total  uncertainty  is  defined  as  the  uncertainty  regarding  the  true 
value  for  an  individual,  selected  at  random  from  the  population.  This  total 
uncertainty  may  be  evaluated  by  the  following  sort  of  two-step  Monte  Carlo 
approach.  (In  pract<ce  th“  'c  /-.-•n  ''fron  he  telescoped  into  one.) 

First,  pick  a  value  at  random  to  represent  the  population  mean,  from  a 
distribution  with  mean  equal  to  the  preferred  value,  and  other  parameters 
(e.g.,  variance)  determined  by  the  variation  that  we  observe  among  studies 
(reflecting  source  M  uncertainty) .  Second,  randomly  select  a  value  to 
represent  the  parameter  value  for  an  Individual  in  the  population,  from  a 
distribution  with  mean  equal  to  the  simulated  population  mean  from  the  first 
step,  and  with  other  parameters  determined  by  the  variation  observed  among 
individuals  within  studies  (reflecting  source  I  uncertainty).  When  the 
complete  process  is  performed  an  indefinite  number  of  times,  performing  each 
step  precisely  once  for  each  final  value  generated,  the  distribution  of  the 
final  values  (from  the  second  step)  describes  the  distribution  of  "total 
uncertainty. " 

Note  that  of  the  two  sources  of  uncertainty,  source  M  is  relatively  more 
subjective  in  nature,  ideally  depending  on  the  scientific  process  of  reviewing 
the  relative  reliability  of  each  study,  and  weighting  studies  appropriately  in 
their  contribution  to  the  preferred  value.  (Uncertainty  may  be  reduced  if  the 
more  extreme  values  are  determined  by  Independent  arguments  to  be  based  on 
relatively  unreliable  procedures  and  are  therefore  given  less  weight.)  In 
view  of  this  process,  source  M  may  be  regarded  as  more  legitimately  the 
subject  of  expert  opinion  than  source  I. 

The  two  basic  sources  of  uncertainty  may  be  subdivided  further  for 
particular  parameters,  depending  upon  the  method  by  which  parameter  estimates 
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were  determined.  For  example,  if  one  assumes  for  humans  a  value  measured 
directly  for  animals,  then  a  component  of  source  M  uncertainty  is  due  to 
variation  among  the  values  typical  for  different  species.  Such  considerations 
can  be  incorporated  into  the  estimation  of  uncertainty  distributions. 

In  many  instances,  it  is  convenient  to  identify  and  estimate  source- 
specific  uncertainty.  Uncertainty  factors  corresponding  to  the  sources  can  be 
denoted  either  as  UFj  or  UF^,  as  appropriate.  These  source-specific 
uncertainty  factors  would  represent  the  residual  uncertainty,  if  uncertainty 
from  other  sources  could  be  eliminated.  Presented  below  are  methods  for 
combining  source  -  specif ic  uncertainty  factors  into  a  total  uncertainty  factor. 

In  the  context  of  sources  of  uncertainty,  it  is  valuable  to  distinguish 
between  reducible  and  irreducible  sources  of  uncertainty.  To  a  large  extent, 
uncertainty  about  the  population  mean  (source  M)  is  reducible  uncertainty,  in 
principle.  Mechanisms  for  reducing  this  uncertainty  include  scrutiny  of 
methods  used  in  particular  studies,  possibly  leading  to  identification  of 
extreme  estimates  that  are  based  on  faulty  procedures.  Also,  reduction  of 
uncertainty  may  be  achieved  by  averaging  estimates  that  are  considered 
comparable  in  reliability.  The  studies  that  are  properly  conducted  and 
equally  reliable  should  all  be  able  to  contribute  to  the  estimation  of  a 
population  average. 

In  contrast.  Individual  variation  as  a  source  of  uncertainty  is  more  or 
less  irreducible.  This  does  depend  on  the  context  of  application  however.  If 
the  goal  of  an  assessment  is  to  describe  the  distribution  of  risk  levels  among 
individuals  in  a  specified  population,  then  some  consideration  of  individual 
variation  in  parameter  values  is  appropriate  and  should  be  Incorporated  into 
the  analysis.  In  that  case  the  object  is  not  to  reduce  the  source  I 
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uncertainty,  but  to  incorporate  a  reliable  estimate  of  that  uncertainty.  On 
the  other  hand,  if  an  assessment  is  desired  for  a  particular  individual,  then 
individual  variation  is  a  relevant  source  of  uncertainty  when  attributes  of 
the  individual  are  unknown.  (In  this  case  we  equate  uncertainty  regarding  the 
individual's  parameter  values  to  the  frequency  distribution  of  values  observed 
for  the  population.)  Obviously  the  uncertainty  may  be  reduced  if  relevant 
characteristics  of  the  individual,  such  as  sex,  weight,  or  breathing  rate,  are 
measured  and  taken  into  account. 

Methods  for  Estimation  of  Uncertainty  Factors.  The  discussion  above  has 
indicated  sources  of  uncertainty  and  the  manner  in  which  uncertainty  about  the 
values  of  a  parameter  can  be  modeled,  i.e.,  either  by  a  lognormal  distribution 
about  a  preferred  value  or  by  a  Dirichlet  distribution  with  expected  values 
for  proportions  set  equal  to  the  preferred  values  of  those  proportions.  What 
remains  is  to  specify  the  methods  that  can  be  used  to  estimate  the  other 
parameters  of  those  distributions,  the  parameters  that  actually  correspond  to 
uncertainty  in  the  sense  of  defining  the  "spread"  of  the  distribution. 

In  the  case  of  the  lognorraally  distributed  uncertainties,  several 
statistics  that  are  computed  directly  from  available  data  are  relevant  to  the 
estimation  of  uncei tainty  factors  (Table  VI-1-17).  Recall  that  for  a 
lognormal  distribution,  the  uncertainty  factor  (UF)  that  is  desired  is  one 
such  that  the  interval  (X/UF,  X*UF)  contains  approximately  95%  of  the  total 
probability  of  the  distribution,  i.e.,  that  specifies  with  about  95%  certainty 
the  interval  that  contains  the  true  value  of  the  parameter,  where  X  is  the 
preferred  value  of  the  parameter.  Three  cases  can  be  identified  that  may  lead 
to  estimates  of  UF  based  on  different  statistics  shown  in  Table  Vl-1-17;  all 
are  Intended  to  yield  UF  estimates  with  the  desired  "95%"  property. 
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In  the  first  case,  suppose  that  data  on  the  variation  of  a  particular 
parameter  is  given  only  in  terms  of  a  plausible  range,  perhaps  specified  by  a 
relevant  expert.  Assume  that  this  range  represents  a  "95%"  Interval  in  the 
sense  given  above.  Then  it  is  natural  to  equate  the  uncertainty  factor  to  the 
statistic  SI  -  (UB/LB)°  where  LB  is  the  lower  bound  of  the  Interval  and  UB 
is  the  upper  bound.  This  representation  can  be  derived  by  equating  LB  to 
GM/UF,  and  UB  to  GM*UF,  where  GM  is  the  geometric  mean  of  the  range. 

The  second  case  arises  when  individual  measurements  (or  relevant  summary 
statistics  of  those  measurements)  are  available  and  it  is  assumed  that  the 
measurements  are  independent  and  identically  distributed  (iid).  The  relevant 
summary  statistics  are  estimates  of  the  log-scale  standard  deviation  (SDL); 
the  primary  ones  considered  are  SSDL  and  CV  (cf.  Table  VI-1-17).  [As 
indicated  in  Table  VI-1-17,  SSDL  is  the  preferred  estimate  of  SDL  but  it  is 
frequently  the  case  that  individual  measurements  are  not  provided,  so  that 
SSDL  can  not  be  computed,  and  that  SSDL  Itself  is  also  not  reported.  Use  of 
CV,  the  coefficient  of  variation,  is  then  sufficient.  Its  approximation  of 
SDL  is  based  on  a  first-order  Taylor's  series  expansion.  Our  experience 
indicates  that  the  approximation  is  accurate  to  one  to  three  significant 
digits.]  Recall  that,  for  a  normal  distribution,  approximately  95%  of  the 
probability  is  contained  within  two  standard  deviations  of  the  mean.  Thus 
for  the  log-normally  distributed  uncertainty  distribution,  the  Interval 
(X/exp(2*SDL) ,  X*exp(2*SDL) )  contains  approximately  95%  of  the  probability. 
When  the  uncertainty  in  question  is  irreducible,  or  when  no  steps  are  taken  to 
reduce  a  reducible  uncertainty,  these  considerations  imply  a  choice  for  UF  of 
S2  (Table  VI-1-17).  When  a  reducible  uncertainty  has  been  reduced  by 
averaging  the  lid  observations  that  underlie  these  statistics,  so  that  the 
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preferred  value  Is  the  geometric  mean  of  the  observations,  then  the  statistic 
S3  (Table  VI-1-17)  appropriately  reflects  the  reduction  in  uncertainty 
attained  by  the  averaging. 

The  third  case  involves  obseir'atlons  that  are  likely  to  deviate  strongly 
from  the  iid  assumption,  with  the  nature  of  the  deviation  not  fully  understood 
or  specified.  This  case  is  obviously  problematic.  This  type  of  situation  may 
hold  frequently  in  evaluating  source  M:  the  available  studies  may  have  used 
different  methods  to  estimate  the  value  of  a  parameter  and  particular  studies 
may  have  associated  with  them  measurement  biases  or  levels  of  measurement 
error  which  are  more  or  less  specific  to  the  study  and  are  of  unknown 
magnitude.  If  some  measure  of  relative  reliability  can  be  attached  to  each 
available  measurement,  and  if  the  relationships  between  the  different 
measurements  can  also  be  fully  described,  then  in  principle  it  may  be  possible 
to  develop  a  formal  assessment  of  uncertainty  using  a  probability  model.  In 
the  absence  of  such  a  complete  analysis  there  is  no  guarantee  that  any  formula 
will  produce  a  better  UF  than  a  guess  from  an  experienced  researcher.  In  some 
cases  one  may  have  to  apply  methods  based  on  iid  assumptions,  leading  to 
evaluations  of  uncertainty  that  may  be  relatively  questionable.  These  cases 
should  be  documented  when  they  occur. 

It  was  argued  above  that  when  the  technique  used  to  reduce  a  reducible 
uncertainty  is  averaging  of  the  relevant  iid  observations,  the  statistic  S3  is 
an  appropriate  estimate  of  the  uncertainty  factor  for  the  preferred  value  that 
is  the  geometric  mean  of  the  observations.  However,  not  all  of  the  preferred 
values  are  geometric  means  of  observations  and,  moreover,  observations  from 
different  studies  (pertinent  to  the  estimation  of  source  M  uncertainty)  are 
frequently  not  lid.  A  development  of  estimates  of  uncertainty  for  other 
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uncertainty- reduction  mechanisms  would  require  a  formal  model  of  such 
mechanisms  that  have  operated  in  producing  the  preferred  values,  including, 
perhaps,  the  impressions  of  experts  regarding  the  relative  merits  of  different 
studies.  This  is  not  a  trivial  undertaking.  However,  it  may  be  the  case  that 
an  adjustment  similar  to  that  seen  in  statistic  S3  (i.e.,  division  by  the 
number  of  measurements)  may  be  considered,  replacing  the  number  of 
measurements  by  some  number  not  larger  than  the  number  of  study- specific 
estimates  considered  in  arriving  at  the  preferred  value.  Again,  source  M 
appears  to  be  more  legitimately  evaluated  on  the  basis  of  expert  opinion  than 
source  I . 

The  statistics  such  as  those  in  Table  VI-1-17  that  are  computed  to 
estimate  source  M  uncertainty  are  all  "impure"  estimates.  The  value  reported 
in  each  study  is  affected  to  some  degree  by  the  variation  of  Individuals 
within  studies,  depending  on  the  number  of  individuals  measured. 

Consequently,  to  an  extent,  individual  variation  is  included  twice  in  such  an 
evaluation  of  total  uncertainty.  Similar  redundancy  occurs  in  evaluating 
various  other  sources  of  uncertainty.  In  principle  a  correction  is  possible 
based  on  a  components -of -variance  approach,  but  such  an  approach  has  not  been 
investigated  here. 

For  those  parameters  that  are  expressed  as  proportions  and  whose 
uncertainty  is  modeled  by  a  Dirlchlet  distribution,  the  uncertainty  parameter 
that  must  be  estimated  is  THETA.  For  proportions  described  by  the  Dirichlet 
distribution,  the  variance  of  an  observed  proportion  f(i)  is  given  by 

(3)  V  -  p(i)*(l-p(l))/(THETA  1), 
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where,  as  discussed  above,  p(i)  is  the  expected  value  of  the  proportion. 
Therefore  a  "quick  and  dirty"  estimate  of  THETA  can  be  obtained  by  equating 
the  variance  of  a  sample  of  measurements  of  a  proportion  to  V,  substituting 
the  mean  proportion  for  p(i).  For  example.  Caster  et  al .  (1956)  reported  that 
the  mean  liver  volume  of  rats,  as  a  proportion  of  total  body  volume,  was 
0.0477,  with  a  variance  of  2.8*10'^.  Thus,  the  approximation  of  THETA  can  be 
determined  from 

(4)  2.8*10'5  -  0.0477*(1  -  0 . 0477)/(THETA  +  1), 


so  that  THETA  is  estimated  to  be  1621.  Other  sample  variances  may  suggest  a 
different  value  for  THETA,  in  which  case  some  sort  of  average  of  the  estimates 
may  be  used.  Figure  VI-1-1  shows  the  Implied  probability  density  function 
for  liver  as  a  proportion  of  body  volume  given  by  the  equation  above.  (Each 
of  the  f(i)  will  follow  a  beta  distribution  with  parameters  p(i)  and  THETA.) 

In  the  discussion  above,  source-specific  uncertainties  (corresponding  to 
source  M  or  source  I)  were  discussed.  The  methods  for  estimating  uncertainty 
factors  just  presented  were  also  specific  to  one  source  of  uncertainty. 

Methods  for  combining  the  source-specific  uncertainty  factor  estimates  to 
obtain  a  total  uncertainty  factor  are  discussed  here. 

In  practice,  when  the  distributions  are  simple,  the  distributions 
associated  with  the  two  sources  can  be  combined  into  one.  Where  the  source-M 
distribution  is  normal  with  mean  Ml  and  variance  VI,  and  the  source-I 
distribution  is  normal  with  mean  0  and  variance  V2 ,  then  the  distribution 
generated  by  the  two-step  algorithm  is  normal  with  mean  Ml  and  variance  V1+V2 . 
This  distribution  approximately  describes  total  uncertainty  for  parameters 
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whose  uncertainty  distribution  is  assumed  to  be  log-normal,  provided  that  the 
means  and  variances  refer  to  the  log  scale. 

In  the  case  of  log-normally  distributed  uncertainty,  if  uncertainty  is 
evaluated  for  K  independent  sources  (e.g.  K-2  for  source  M  and  source  I)  then 
we  may  estimate  a  UF  specific  to  each  source  (UF^  for  source  i,  i  -  1,...,K.). 
Corresponding  to  UF^  is  a  log-scale  standard  deviation  SDLj  -  ln(UFi)/2  for 
the  assumed  lognoirmal  probability  distribution.  If  the  sources  of  uncertainty 
are  assumed  to  be  Independent,  then  the  appropriate  value  of  SDL  representing 
our  total  uncertainty  (SDLj)  is 

(5)  SDLi  -  (El  SDLi^  ]°  ^ 


Consequently  an  appropriate  UF  for  total  uncertainty  is 


(>  UF,  -  exp({Ei(lnUFi)2  )°-5]. 


Thus,  for  example,  any  UFi  equal  to  unity  (l.e.  any  SDLi  equal  to  zero)  makes 
no  contribution  to  the  overall  UF. 

These  arguments  apply  to  uncertainty  associated  with  use  of  direct 
measurements  of  a  quantity  of  Interest.  In  some  cases,  for  example  when  a 
measurement  of  a  parameter  for  a  rodent  is  used  for  humans,  special  treatment 
may  be  required.  In  any  case,  uncertainty  can  be  evaluated  by  constructing  an 
algorithm  such  as  the  "two-step"  algorithm  above,  where  the  first  step  is  to 
pick  a  population  mean  value  from  a  distribution  with  mean  equal  to  the 
preferred  value  and  variance  reflecting  the  source  M  uncertainty,  and  the 
second  step  is  to  pick  a  value  for  the  parameter  from  a  distribution  having  a 
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mean  equal  to  the  value  s  .^pled  in  step  one  and  a  variance  reflective  of 
individual  vaT-iatlon.  For  the  example  just  cited,  che  variance  for  the  first 
step  must  include  consideration  of  the  differences  in  the  parameter  observed 
for  different  species. 

It  is  somecines  the  case  that  a  PEPK  model  parr^meter  is  not  measured 
directly  but  is  a  function  of  other  variables  that  are  themselves  su’  ject  to 
uncertainty.  Let  f(h)  be  a  function  of  the  vector  of  variables 
h'  -  (hi,...,hK),  each  subject  to  uncertainty.  For  example,  a  tlssue/blood 
partition  coefficient  (PC)  for  a  given  compound  is  generally  estimated  by  the 
ratio  of  an  estimate  of  the  tissue/gas  partition  coefficient  (here  denoted  h^) 
to  an  estimate  of  the  blood/gas  partition  coefficient  (here  denoted  h2) , 
i.e.,  PC  -  t(h)  -  hi/h2.  It  is  nece.ssary  to  evaluate  uncertainty  in  the 
function  f  implied  by  uncertainty  in  its  arguments  hi  and  h2. 

The  simplest  case  is  when  the  function  is  linear.  For  example,  in  the 
example  of  partition  coefficients,  it  is  convenient  to  deal  with  uncertainty 
in  the  natural  log  of  the  tissue/blood  PC,  since  we  will  assume  that  the  PC's 
are  normally  distributed  in  the  log  scale,  with  specif it d  variance  SDL2. 
Therefore  we  evaluate  uncertainty  for  the  function 

(7)  In  PC  -  In  hi  -  In  h2. 


A  further  simpllcation  is  to  consider  the  uncertainties  in  the  arguments  to  be 
independent.  In  this  case  the  variance  of  a  linear  function  is  simply  the 
same  linear  function  of  the  variances  of  the  arguments,  but  with  the 
coefficients  squared  (again,  assuming  that  the  arguments  are  distributed 
independently.;  If  we  assume  independent  lognormal  distributions  for  hi  and 
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hj,  with  log-scale  standard  deviations  SDLj  and  SDL2,  respectively,  then  the 
log- scale  standard  deviation  for  the  function  f  is 


|~  17)  SDIf  -  [  (SDLi^  SDLa^)  jO.s 


This  is  essentially  formula  (5).  The  app'-opriate  UF  is  simply  exp(2*SDLf)  . 

For  'linear  functions  of  uncertain  arguments,  simple  Taylor-series 
approximations  .re  available;  also,  non- >idependent  uncertainties  in  the 
arguments  can  bt  handled  in  a  simple  manner,  provided  that  the  dependencies 
are  well  described  by  cor  Intlons  or  covariances. 

Uncertainty  Evaluation  for  PERC  Partition  Coefficients.  The  first 
example  of  hew  uncertainty  from  various  sources,  or  arising  from  various 
assumptions,  can  be  combined  into  an  assessment  of  total  uncertainty  for  PEIRC 
PBPK  model  parameters  was  applied  to  the  partition  coefficients.  A  partition 
coefficient  (PC)  is  the  ratio  of  concentrations  in  two  compartments  at 
equilibrium,  and  so  is  a  dimensionless  quantity  that  is  logica'’ly  bounded 
below  by  zero  but  not  logically  bounded  above.  In  ac^'ordance  with  the 
discussions  above,  the  uncertainty  was  modeled  using  a  lognormal  distribution 
with  log-scale  mean  equal  to  the  log  of  the  preferred  value. 

The  examples  of  uncertainty  estimation  providec  here  illustrate  four 
different  scenarios.  First  were  those  PC's  that  were  measured  directly,  i.e. 
the  blood/a) r  PC's,  and  which  may  have  been  measured  by  several  investigators. 
Second,  for  the  rat  tlssue/blood  PC's,  the  estimates  were  derived  as  functions 
(ratios;  of  tissue/air  and  blood/air  PC's  and  thus  their  uncertainty  was 
evaluated  u.sing  equation  (7).  Third,  human  tlssue/air  PC's  were  not  measured 
directly  so  that  rodent  values  had  to  be  attributed  to  humans,  lllu..trstlng 
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the  Importance  of  considering  species- to-species  variation  In  the  estimation. 
Fourth  were  the  PC's  such  as  rapidly  perfused/air  that  were  not  measured  in 
any  species  but  were  equated  with  specific  other  PC's,  in  this  case  the 
liver/atr  PC.  The  derivation  of  uncertainty  estimates  is  illustrated  for  each 
case  . 

For  all  of  these  "estimation  scenarios"  source  M  uncertainty  was 
visualized  as  shown  in  Figure  VI-1-2.  For  each  tissue  group,  there  is  a 
(hypothetical)  super-distribution  of  values  of  PC's  for  that  group  over 
mammalian  species.  The  species  of  greatest  interest  in  the  present  context 
were  rats,  mice,  and  humans.  Each  species  has  for  each  tissue  group  its  own 
distribution  of  values  of  PC's  and  an  associated  true  mean  value.  This 
distribution  is  manifested  in  the  measurements  that  are  obtained  from 
individual  animals  and  the  mean  of  those  measurements  from  each  study.  Source 
I  uncertainty  can  be  visualized  in  terms  of  additional  branches  added  to  the 
top  of  Figure  VI-1-2,  branches  emanating  from  the  means  of  the  measured  values 
to  the  measurements  taken  from  individual  animals. 

The  representation  given  by  Figure  VI-1-2  has  associated  with  it  certain 
assumptions.  The  first  is  that  a  tissue  group  in  one  species  is  more  similar 
to  that  tissue  group  in  another  species  than  that  tissue  group  is  to  other 
tissue  groups  in  the  same  species.  This  appears  to  be  supported  by  the 
observations  of  Fiserova-Bergerova  (1983)  on  species  differences  in  partition 
coefficients.  She  stated  that  species  differences  in  PC's  appear  to  be 
haphazard  and  most  likely  without  biological  significance.  If  this  is 
realistic,  then  perhaps  it  is  appropriate  to  treat  species- to-specles 
variation  in  PC's  as  independent  for  different  compartments.  This  is  the 
second  assumption  underlying  the  representation  given  in  Figure  VI-1-2. 
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As  Indicated  In  the  figure,  the  estimation  of  the  true  mean  of  a  species- 
specific  distribution  requires,  at  least,  the  estimation  of  that  mean  from  the 
measured  mean  values.  This  has  associated  with  it  uncertainty  represented  in 
the  figure  by  SMl(x,i),  a  log-scale  standard  deviation  for  species  x  and 
tissue  group  1  that  reflects  the  variation  one  would  expect  to  see  in  ""ean 
values  given  the  distribution  around  the  true  mean  PC  for  that  species  and 
tissue  group.  Other  components  of  source  M  uncertainty  are  displayed  in  the 
figure  as  SM2(i)  and  SM3(i,J),  log-scale  standard  deviations  representing 
variation  across  species  of  mean  values  of  PC^  and  variability  in  the 
relationship  between  PC^  and  PCj ,  respectively.  Specific  estimates  for  each 
of  the  standard  deviations  are  discussed  below,  but  it  is  perhaps  appropriate 
to  note  here  that,  given  those  estimates,  the  overall  source  M  uncertainty 
(expressed  in  terms  of  a  lognormal  variance)  can  be 

determined  by  summing  the  squared  estimates  of  the  standard  deviations  that 
are  associated  with  each  branch  that  must  be  traversed  to  get  from  measured 
means  to  the  true  mean  desired. 

Estimates  of  uncertainties  for  the  directly  measured  PC's,  the  blood/alr 
coefficients,  were  considered  first.  In  terms  of  Figure  VI-1-2,  the 
uncertainty  Involved  is  associated  with  the  single  branch  from  the  measured 
values  to  the  true  mean.  Tables  VI-1-18  and  VI-1-19  summarize  blood/air 
partition  coefficient  estimates  that  have  been  obtained  for  various  chemicals, 
including  PERC.  Source  M  uncertainty  in  this  case  is  identified  solely  with 
variation  across  studies.  Thus,  as  indicated  in  Table  VI-1-19,  if  no 
uncertainty  reduction  is  attempted,  the  value  0.32  Is  the  estimate  of  SMl(h,b) 
and  the  statistic  S2 ,  equal  to  1.9  for  PERC  in  humans,  is  the  estimate  for 
UF^.  However,  if  the  geometric  mean  of  the  four  PERC  human  blood/alr  PC 
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estimates  is  used  as  the  preferred  value,  and  If  the  measurements  from 
different  studies  can  be  assumed  to  be  iid  (a  questionable  assumption),  then 
the  estimate  of  SMl(h,b)  (tor  the  mean  of  the  measured  values)  Is  0. 32/(4)®  * 
-0.16  and  UFm  may  be  equated  to  S3,  which  In  this  case  Is  equal  to  1.4. 

No  direct  information  is  available  on  the  individual  variation  In  the 
human  PERC  blood/air  PC  for  estimation  of  UFj .  However,  Table  VI-1-20 
displays  some  information  on  individual  variation  for  other  volatile 
lipophilic  compounds  that  may  be  used  in  the  absence  of  more  pertinent  data. 
That  table  indicates  that  inter- individual  variation  of  the  blood/alr  PC  is  at 
a  maximum  for  trichloroethylene  in  humans.  Thus,  let  us  assign  the  SOLE 
estimate  for  that  case  to  the  SOLE  that  will  be  used  to  characterize  inter¬ 
individual  variability  for  the  PERC  blood/air  PC  in  humans.  Using  statistic 
S2  to  estimate  UFj ,  one  obtains  the  value  1.8  (-exp(2*0. 3)) .  Finally, 
combining  UF^  and  UFj  via  equation  (6),  the  estimate  for  UF^  is 

(8)  UFx  -  exp[  {ln(1.4)2  +  in(1.8)^)0  *]  -  2.0. 


If  the  same  considerations  were  applied  to  rats  and  mice,  then  the  first 
immediate  problem  would  be  that  there  are  no  data  regarding  study -to -study 
differences  in  blood/air  PC's  for  these  two  species.  The  data  in  Tables 
VI- 1-18  and  VI -1-19  were  used  to  address  this  issue.  We  tentatively  based 
evaluations  for  various  tlssue/air  PC's  on  data  for  blood/air  PC's,  assuming 
similar  levels  of  relative  variation  between  studies  for  different  tissue/air 
PC's.  Also,  we  know  of  no  reason  to  suspect  different  levels  of  relative 
variation  between  studies  for  different  blood/air  PC's,  provided  that  the 
values  of  the  PC's  themselves  are  not  too  different  in  magnitude.  (Based  on 
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preliminary  investigations,  very  small  PC's  were  expected  to  be  accompanied  by 
somewhat  larger  SDL's,  representing  measurement  error;  consequently,  we 
ignored  the  large  values  of  S2  in  Table  VI-1-19  which  correspond  tc  small 
values  of  the  GM,  e.g.,  9.8  for  carbon  tetrachloride  and  3.4  for  ethane.) 

Based  on  the  S2  column  in  Table  VI-1-19,  it  appeared  that  a  UF  of  2  was  a 
reasonable  (or  slightly  inflated)  estimate  for  source  M,  for  many  tissue/air 
PC's,  if  uncertainty-reducing  procedures  are  not  performed  or  taken  into 
account.  In  particular,  a  of  2.0  for  rat  or  mouse  blood/air  PC's  was 
assumed.  Expressed  in  terms  of  the  standard  deviations  displayed  in  Figure 
VI-1-2,  that  estimate  of  UF^  corresponds  to  SMl(r,b)  -  SMl(m,b)  -  0.347. 

Similarly,  the  data  in  Table  VI-1-20  were  used  in  the  absence  of  data  on 
inter- individual  variation  in  blood/air  PC  values  in  rodents.  The  maximal 
SOLE  for  the  blood/alr  PC  from  that  table  was  0.3,  the  same  value  used  in  the 
derivation  of  the  UFj  for  human  blood/alr  PC.  Thus,  the  same  UFj  value  was 
estimated  for  the  rodent  blood/air  PC's  as  for  the  human  PC,  1.8.  Combining 
UF^  and  UFi  via  equation  (6) ,  one  obtains  a  total  uncertainty  factor  for 
rodent  blood/air  PC's  as  follows: 

(9)  UFj  -  exp[  {ln(2.0)2  ln(l .  8)2}0-5]  _  2.5. 


The  second  case  involved  those  rodent  PC's  that  are  ratios  of  tissue/air 
and  blood/air  PC's.  The  derivation  of  an  uncertainty  factor  in  these 
instances  involved  the  use  of  equation  (7).  An  example  for  liver/blood  PC  in 
rats  is  provided. 

Once  again,  no  study- to -study  variability  data  were  available  for  rats 
for  this  PC.  As  argued  above,  the  data  in  Table  VI- 1-19  suggested  that  UF^  - 
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2.0  was  not  unreasonable,  In  this  case  both  for  liver/air  and  blood/air  PC's. 
The  application  of  equation  (7)  to  the  two  (Identical)  UF^  values  yields  a 
velue  of  2.7  for  for  rat  liver/blooH  pc. 

Table  VI -1-20  presents  data  on  the  individual  variation  in  the  value  of 
liver/air  partition  coefficients.  No  data  on  rats  were  available,  but  the 
maximum  value  for  SOLE  corresponding  to  the  liver/air  PC  was  1.0.  The  use  of 
the  maximum  SOLE  is  recommended  for  all  of  those  cases  in  which  the  partlcul cr 
species  of  interest  is  not  measured  directly.  In  this  instance,  the  maximum 
value  was  from  the  study  of  Webb  and  Weaver  (1981)  using  horses.  The  inter- 
individual  variation  for  those  horses  is  likely  to  be  greater  than  that  for 
carefully  bred  laboratory  rats.  The  SOLE  estimate  of  1.0  corresponds  to  a  UFj 
estimate  for  the  liver/air  PC  of 

(10)  exp(2*1.0)  -  7.4. 


Using  equation  (7)  to  combine  Source  1  uncertainties  for  the  liver/air  and 
blood/air  PC's,  one  obtains  UFj  -  8.1  for  the  rat  liver/blood  PC.  Combining 
UFf,  and  UFj  using  equation  (6),  the  resulting  UFj  was 

(11)  UFt  -  exp[ {ln(2.7)2  +  ln(8.1)2}0  =]  -  10.1. 

I _ 

A  completely  analogous  situation  applied  to  the  derivation  of  rat  fat/blood 
and  slowly  perfused/blood  PC  uncertainty  factors  (if  the  slowly  perfused 
tissues  are  equated  with  muscle). 

The  third  case  for  PC  uncertainty  estimation  involved  those  PC's  that 


were  not  measured  in  the  species  for  which  estimates  are  desired.  Human 


tissue/alr  PC's,  and  often  mouse  tissue/air  PC's,  have  not  been  measured 
directly  for  tissues  other  than  the  blood.  The  use  of  the  PC's  estimated  from 
rat  tissues  to  approximate  the  human  or  mouse  values  added  uncertainty  to  the 
human  or  mouse  estimates  that  was  related  to  the  variation  in  PC's  seen  across 
species.  In  terms  of  Figure  VI -1-2,  this  uncertainty  Involved  the  branches 
connecting  the  measured  species -specific  mean  values  and  the  two  branches 
connecting  the  tissue  group  to  the  specific  species  of  Interest,  the  one  for 
which  measurements  were  available  and  the  one  for  which  an  estimate  of  the 
population  mean  was  desired.  An  example  using  the  human  liver/blood  PC  is 
presented. 

The  liver/air  PC  value  measured  for  rats  was  to  be  attributed  to  humans. 
It  was  shown  above  that  the  source  M  uncertainty  associated  with  the  rat 
estimate  itself  yields  -  2.0  (equivalently,  SMl(r,l)  -  0.347).  This 
uncertainty  had  to  be  combined  with  the  uncertainty  due  solely  to  the  need  to 
extrapolate  the  estimate  across  species.  That  is,  SM2(1)  had  to  be  estimated 
also.  The  resulting  equation  for  combined  uncertainty,  expressed  in  terms  of 
the  standard  deviation  is  as  follows: 

(12)  SOLE  -  (2*SM2(1)2  +  SMl(r , l)2)0-5, 


where  the  branches  corresponding  to  the  need  to  go  from  the  true  mean  rat 
liver/air  PC's  to  the  true  mean  human  liver/alr  PC  contributed  the  associated 
standard  deviation  term,  SM2(1),  twice. 

Table  VI-1-21  gives  data  from  Flserova-Bergerova  (1983)  on  variation  in 
tissue/alr  PC's  among  four  species  (human,  monkey,  dog,  and  rat)  for 
isoflurane  and  methylene  chloride.  Some  statistics  based  on  this  data  are 
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reported  in  Table  VI-1-22.  Since  the  estimate  used  by  Reitz  and  Nolan  (1986; 
the  source  of  preferred  values  considered  here)  Involved  attributing  rodent 
estimates  to  humans,  the  most  relevant  data  from  Table  VI-1-21  appeared  to  be 
that  for  rats  and  humans.  Therefore,  Included  in  Table  VI-1-22  were  estimates 
of  the  error  to  be  expected  in  extrapolations  of  this  type  ("human/rat 
error").  It  was  reassuring  that  the  uncertainty  factors  estimated  by 
statistic  S2  were  at  least  as  large  as  the  human/rat  error  term.  Some 
additional  conservatism  may  be  warranted  since  the  ild  assumption  Involved  in 
equating  a  UF  to  S2  may  be  doubtful  in  this  case. 

Based  on  the  values  of  S2  shown  in  Table  VI-1-22,  use  of  a  UF  of  1.8 
appeared  reasonable  or  slightly  conservative  for  most  tissues.  (An  exception 
was  muscle,  for  which  S2  was  2.3  for  isoflurane  and  2.0  for  methylene 
chloride.)  Thus,  an  SM2(1)  based  on  a  UF  of  1.8  for  the  liver  compartment  was 
used;  in  this  case  SM2(1)  was  estimated  to  be  ln(1.8)/2  -  0.29.  Combining 
this  SM2(1)  with  SMl(r,l)  (-0.34)  via  equation  (12)  yielded  a  standard 
deviation  of  0.54  or  an  uncertainty  factor  of  2.92. 

The  calculations  just  presented  were  for  llver/air  partitioning.  To  get 
the  liver/blood  PC  for  humans,  equation  (7)  was  used  to  combine  UF^'s  for 
liver/alr  and  blood/air.  Recall  that  for  humans  the  blood/air  PC  UF^  was 
equal  to  1.4.  Thus  equation  (7)  yielded 

(13)  UF„  -  exp(  {ln(2.9)2  -H  ln(1.4)*}'>-5)  -  3.1 


for  the  human  liver/blood  source  M  uncertainty  factor. 

For  source  I,  the  data  on  PERC  PC's  in  humans  were  as  scarce  as  those 
for  PERC  PC's  in  rats.  Thus  the  same  procedure  as  described  above  for  the  rat 
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UFj  calculations  was  followed.  The  maximum  Inter- Individual  variation  was 
used  (cf.  Table  VI-1-20).  The  SOLE  of  1.0  was  consistent  with  UFj  -  8.1. 

Total  uncertainty  was  reflected  In  the  value  UFj  -  10.8  derived  using  equation 
(6)  with  the  indicated  values  of  UF^^  and  UFj . 

The  uncertainties  associated  with  the  human  fat/blood  and  slowly 
perfused/blood  PC's  were  estimated  in  an  analogous  manner. 

The  final  situation  that  arises  in  the  estimation  of  partition 
coefficients  is  Illustrated  by  the  rapidly  perfused/blood  PC  (also  referred  to 
as  the  rapid/blood  PC)  in  rats  and  humans.  In  this  case,  no  direct 
measurements  have  been  taken  for  this  tissue  group;  instead,  the  value  of  its 
PC  was  inferred  from  measurements  on  other  tissue,  l.e.,  the  liver.  Once 
again,  this  is  an  uncertainty  associated  with  source  M,  i.e.,  with  the 
estimation  of  a  population  mean  value. 

The  estimation  of  a  PC  for  one  group  of  tissues  using  PC  measurements 
from  another  tissue  is  problematic.  Although  the  liver  may  be  representative 
of  rapidly  perfused  tissue  in  many  respects,  it  is  not  at  all  certain  that  a 
PC  for  liver  is  representative  of  a  PC  characterizing  the  entire  tissue  group. 
Indeed,  as  the  data  in  Table  VI-1-20  show,  other  rapidly  perfused  tissues  such 
as  the  brain,  kidney,  or  lung  may  have  PC's  that  differ  substantially  from 
that  of  the  liver.  If  one  represents  the  estimation  of  the  rapid/air  PC  from 
the  liver  PC  by 

(14)  PCj  -  e*PCi, 


where  a  is  a  coefficient  that  must  be  estimated,  then  the  uncertainty  arises 
in  the  estimation  of  that  coefficient.  Reitz  and  Nolan  (1986)  assumed  that  a 
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was  equal  to  1.  The  data  in  Table  VI -1-20  indicate  that  a  value  for  c  of 
around  0.5  to  0.7  may  be  reasonable  for  some  of  the  chemicals  shown  there. 
However,  PERC  itself  is  not  represented  in  that  table,  and  the  variation  seen 
in  that  table  is  substantial. 

Based  solely  on  the  information  available  in  that  table,  an  uncertainty 
factor  for  this  aspect  of  the  estimation  of  the  rapld/air  PC  of  2  appeared 
justified  and  was  used.  If  the  rapid/air  PC  is  set  equal  to  half  of  the 
liver/air  PC,  then  this  uncertainty  factor  has  the  effect  of  making  the 
Itver/air  PC  the  95%  upper  bound  (of  the  uncertainty  distribution)  for  the 
rapid/air  PC. 

This  cross-tissue  extrapolative  uncertainty  is  reflected  in  Figure  VI-1- 
2  by  the  dotted  line  connecting  the  different  tissue  groups.  It  must  be  added 
to  the  other  uncertainties  that  are  pertinent  to  the  estimation  of  the 
rapid/air  PC,  those  that  are  reflected  in  the  SMI  and  SM2  standard  deviation 
terms.  The  formula  that  combines  the  branch- specif ic  variance  estimates 
pertinent  to  this  estimation  problem  is 

(15)  SOLE  -  (SMl(r,l)2  +  SM2(1)2  +  SM2(r)2  +  SM3(1 , r)^]®-’, 


where  SM3(l,r)  is  the  standard  deviation  associated  with  the  extrapolation 
across  tissue  groups  just  discussed.  Both  SM2  values  were  set  equal  to  0.29 
on  the  basis  of  the  data  in  Tables  VI-1-21  and  VI-1-22  (discussed  above).  The 
value  of  SM3(l,r)  was  glvan  by 

(16)  SM3(l,r)  -  ln(2.0)/2  -  0.347, 
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where  the  value  2.0  was  the  uncertainty  factor  postulated  for  the  cross-tissue 
extrapolation.  Thus,  equation  (15)  yielded  an  estimate  of  0.64  for  the  SOLE 
of  interest.  This  corresponds  to  an  uncertainty  factor,  UFf,.  of  3.6  for 
rapid/alr  partitioning.  Combining  this  with  the  source  M  uncertainty  for 
blood/alr  PC's  (2.0  for  rats,  1.4  for  humans)  using  equation  (7),  the  UFj,'s 
for  rapld/blood  PC's  were  4.3  and  3.8  for  rats  and  humans,  respectively. 

Source  I  uncertainty  was  again  estimated  on  the  basis  of  the  data  shown 
in  Table  VI -1-20.  The  brain  and  the  kidney  were  assumed  to  be  representative 
of  rapidly  perfused  tissues  in  general  with  respect  to  inter- individual 
variation  in  PC  values.  The  brain  SOLE  was  0.5;  the  SOLE  for  kidney  was  only 
slightly  smaller,  0.4.  We  assumed  that  the  tissue  group  as  a  whole  Is  not 
more  variable  than  the  most  variable  component,  in  this  case  the  brain.  An 
SOLE  of  0.5  corresponds  to  UFj  -  2.7  for  the  rapid/air  PC.  Using  equation  (7) 
to  combine  this  uncertainty  with  that  for  blood/air  source  I  uncertainty  (UFj 
-  1.8)  the  resulting  UFj  was 

(17)  UFi  -  exp[  {ln(2.7)2  +  ln(l . S)^}'’-’)  -  3.2. 


Finally,  combining  the  source  M  and  source  I  uncertainties  via  equation 
(6)  gave  the  total  uncertainty  factors,  UFj  -  6.5  for  rats  and  UFi  -  5.9  for 
humans . 

This  concludes  the  examples  of  uncertainty  estimation  for  the  partition 
coefficients  of  the  PERC  model.  The  other  coefficients  were  analyzed  in  the 
same  manner  as  the  examples  just  presented.  The  entire  process  is  summarized 
in  Table  VI-1-23. 
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The  entire  set  of  partition  coefficient  estimates  (preferred  values)  and 
their  estimated  uncertainty  factors  are  displayed  in  Table  VI-1-24.  Note  that 
the  preferred  values  for  human  i/blood  PC's  (where  1  is  one  of  the 
compartments  used  in  the  model)  were  based  on  an  average  of  mouse  and  rat 
i/air  PC's.  The  uncertainty  factors  derived  did  not  consider  this  averaging 
and  any  corresponding  reduction  in  uncertainty.  That  is,  as  the  examples 
above  illustrated,  the  human  PC  uncertainty  factors  were  derived  as  if  only 
one  rodent  species  formed  the  basis  of  the  human  coefficients.  This  is 
probably  not  too  bad.  Even  in  the  lid  case,  the  log-scale  standard  deviations 
that  would  correspond  to  the  averages  would  be  the  unaveraged  standard 
deviations  divided  by  the  square  root  of  2  (the  number  of  terms  contributing 
to  the  average).  That  would  be  a  small  change  in  and  of  itself.  Moreover, 
since  the  iid  assumption  is  at  best  questionable  in  this  case  of  two  sets  of 
rodent  PC's,  the  appropriate  adjustment  should  probably  fall  between  1  and  the 
square  root  of  2.  Because  the  exact  value  of  the  appropriate  adjustment 
factor  is  not  known,  and  because  it  will  make  little  difference  in  the 
estimation  of  uncertainty  factors,  the  unadjusted  uncertainty  estimates  are 
presented  and  are  considered  adequate. 

Estimation  of  Uncertainty  for  PERC  Metabolic  Constants.  The  parameters 
vmax  and  km  define  metabolism  for  the  PERC  PBPK  model.  Values  for  these 
parameters  were  reported  by  Reitz  and  Nolan  (1986)  and  Hattis  et  al.  (1986). 
Reitz  and  Nolan  employed  an  optimization  technique  to  estimate  the  values 
appropriate  for  mice,  rats,  and  humans.  They  used  the  data  of  Schumann  et  al. 
(1980)  and  gas  uptake  data  obtained  from  WPAFB  to  derive  two  sets  of  mouse 
values.  The  data  of  Pegg  et  al .  (1979)  were  used  to  optimize  the  rat  values 
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and  the  data  of  Monster  et  al.  (1979)  and  Fernandez  et  al,  (1976)  were  used 
for  h'unan  optimization. 

Hattls  et  al .  did  not  employ  a  formal  optimization  procedure.  They  used 
the  same  data  sets  for  rat  and  mouse  parameter  estimation  as  used  by  Reitz  and 
Nolan.  Hattls  et  al.  also  referred  to  the  results  presented  by  Mltoma  et  al. 
(1985)  for  parameter  estimation  in  the  two  rodent  species.  For  their  human 
metabolic  constant  estimates,  Hattls  et  al.  adjusted  vmax  and  km  based  on 
Ohtsuj 1  et  al .  (1983)  and  Ikeda  et  al.  (1972).  Described  here  are 
considerations  of  uncertainty  related  to  the  metabolic  parameter  estimates. 

First,  some  comment  on  the  use  of  optimization  to  derive  parameter 
estimates  is  warranted.  Typically,  a  subset  of  PBPK  parameters  is  estimated 
by  optimizing  the  fit  of  model  predictions  to  data,  while  the  remaining 
parameters  are  fixed  at  constant  values.  In  general  the  "fixed"  subset  will 
include  parameters  fo'-  which  in  vitro  measurements  are  available;  the  in  vitro 
values  are  assumed  equal  to  in  vivo  values.  The  choice  of  which  parameters 
will  be  in  each  subset  is  not  necessarily  clear-cut.  For  example,  Reitz  and 
Nolan  fixed  values  of  partition  coefficients  for  their  "Mouse  1"  estimates, 
but  optimized  those  values  for  their  "Mouse  2"  estimates;  vmax  and  km  were 
fitted  for  both  Mouse  1  and  Mouse  2  sets. 

Estimates  for  the  fitted  parameters  depend  to  some  degree  on  the  values 
assumed  for  the  fixed  set.  Uncertainty  for  the  fitted  set  will  be  determined 
by  uncertainty  in  the  fixed  set,  as  well  as  by  experimental  error  in  the  data 
to  which  the  model  is  fit,  by  variations  in  the  results  of  experiments  used  in 
optimization,  and  by  uncertainty  regarding  model  structure.  In  principle, 
uncertainty  regarding  the  values  of  the  fixed  parameters  can  be  evaluated  by  a 
Monte  Carlo  approach  in  which  values  are  selected  randomly  for  the  fixed  set, 
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re-fitting  the  model  with  each  selection.  That  approach  would  be  high'/ 
computer- Intensive ,  For  each  selection  of  values  for  the  fixed  parameters  (a 
selection  that  would  have  to  be  repeated  many  times  to  get  an  accurate  picture 
of  the  variation)  an  optimization  would  have  to  be  performed;  the 
optimizations  themselves  require  many  simulation  runs. 

Rather  than  adopt  such  a  ccmp*,ex  procedure,  we  obtained  a  tentative 
evaluation  of  uncertainty  by  treating  the  available  estimates  of  the  fitted 
parameters  as  a  representative  sample  of  the  values  that  were  reasonable. 
Implicitly,  the  vectors  of  fixed  parameters  associated  with  the  various  fitted 
estimates  were  also  tre -."ed  as  a  representative  sample. 

The  units  used  to  report  \Tiiax  and  km  values  in  Hattis  et  al .  (1986)  were 
not  the  same  as  those  used  by  Reitz  ai.d  Nolan  (1986).  The  Hattis  et  al . 
units  were  converted  to  those  used  by  Reitz  and  Nolan,  using  the  standard 
body  weights  assumed  by  Hattis  et  al .  (0.025,  0.25,  and  70  kg  for  mouse,  rat, 
and  hu..ici,  respectively).  Also,  since  Hattis  et  al  .  defined  metabolism  in 
terms  of  liver  concentration  (as  opposed  to  the  use  of  concentration  in  the 
venous  blood  leaving  the  liver,  the  formulation  employed  by  Reitz  and  Nolan) 
their  km  values  were  converted  to  units  comparable  to  those  of  Reitz  and  Nolan 
by  dividing  by  pi,  the  li.ver/blocd  partition  coefficient.  (This  follows  from 
the  fact  th  t .  according  to  the  models,  liver  concentration  equals  hepatic 
portal  venous  concentration  times  pi.)  In  making  the  conversion,  the  "best 
estimates"  of  pi  reported  by  Hattis  et  al.  were  used  (4.73  for  humans,  7.7_ 
for  rodents).  Vmax  was  .scaled,  as  was  done  in  Reitz  and  Nolan,  by  estimating 
vmaxc  such  that 
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(18)  vmax  -  vmaxc  * 


All  uncertainty  estimates  were  estimated  for  vmaxc  and  km  as  defined  by  Reitz 
and  Nolan. 

Tables  VI-1-25  and  VI-1-26  present  the  estimates  reported  in  Reitz  and 
Nolan  (1986)  and  Hattis  et  al .  (1986).  Hattis  et  al .  gave  two  estimates  for 
each  species,  based  on  different  data  sets,  for  both  km  and  vmaxc.  The 
estimates  were  named  according  to  the  sources  of  data  used  tc  estimate  the 
parameters.  Note  that  the  "Mouse  1”  and  human  estimates  of  km  (9.56)  reported 
in  Reitz  and  Nolan  were  actually  extrapolated  from  rats,  so  they  were  not  used 
for  variance  estimation  in  mice  or  humans. 

For  the  different  estimates  of  metabolic  parameters  reported  in  Tables 
VI-1-25  and  VI-1-26,  the  iid  assumption  c'^errly  seemed  to  be  violated. 
Consequently,  methods  based  on  iid  assumptions  were  applied  only  because  of  a 
lack  of  additional  information. 

For  vmaxc,  values  of  the  sample  standard  deviation  of  logs  (SSDL)  are 
0.195  for  humans,  0.919  for  rats,  and  0.721  for  mice.  These  differences  do 
not  necessarily  reflect  real  differences  between  species  In  the  uncertainty 
involved  in  measuring  vmoxc ,  since  SSDL  is  itself  a  statistic  subject  to  some 
measurement  error.  Tentatively  we  based  an  uncertainty  factor  for  vmaxc  on  an 
SSDL  of  0.532,  tie  weighted  average  of  SSDL' s  from  different  species.  Some 
reduction  of  the  uncertainty  may  be  appropriate  if  one  used  the  average  of 
several  estimates  as  the  preferred  value  for  each  species;  however,  proper 
application  of  the  statistic  S3,  which  accomplishes  such  a  reduction  in  a 
formal  manner,  is  based  on  the  lid  assumption.  Some  conservatism  (implying 
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greater  uncertainty)  was  introduced  by  using  S2  rather  than  S3.  Accordingly, 
the  derived  uncertainty  factor  was  2.9  (-exp [ 2*0. 532 ] )  for  mice,  rats,  and 
humans.  (Uncertainty  could  still  be  reduced  when  S2  is  used,  but  by  a 
different  mechanism:  if  the  spread  of  values  that  were  considered  reasonable 
estimates  was  lowered  by  further  study  of  estimation  procedures,  this  would  be 
reflected  in  a  smaller  value  of  S2.) 

Uncertainty  for  km  was  evaluated  in  a  similar  manner,  arriving  at  an 
uncertainty  factor  of  13  (-exp [ 2*1 . 28] )  for  all  three  species. 

Tables  VI-1-25  and  VI-1-26  also  display  the  geometric  means  of  vmaxc  and 
km  for  all  three  species.  In  the  absence  of  a  reliable  analysis  of  the 
relative  merits  of  the  different  estimates,  we  adopted  the  geometric  means  as 
preferred  values  for  the  Monte  Carlo  simulations.  Note  the  differences 
between  the  geometric  means  and  the  estimates  of  Reitz  and  Nolan. 

The  uncertainty  factors  for  vmaxc  and  km  appeared  to  be  representing 
source  M.  We  suspected  that  individual  variation  (source  I)  was  small  in 
proportion  to  the  source  M  uncertainty  represented  by  factors  of  2.9  for  vmaxc 
and  13  for  km.  In  any  case,  quantitative  Information  suitable  for  evaluating 
individual  variation  for  vmaxc  and  km  was  not  available.  Thus,  tentatively, 
we  set  the  UFj's,  representing  total  uncertainty,  equal  to  the  UF(,'s  of  2.9 
for  vmaxc  and  13  for  km. 

The  statistical  procedure  for  evaluating  individual  variation  in 
metabolic  parameters  estimated  by  optimization  depends  on  whether  the 
available  data  consist  of  multiple  measurements  at  different  times  for  the 
same  individuals,  or  (e.g..  Young  and  Wagner,  1979)  different  animals  measured 
at  each  time.  In  the  first  case,  the  optimization  may  be  performed  separately 
for  each  subject,  and  the  variation  of  subject-specific  estimates  summarized. 
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In  the  second  case  the  model  may  be  fitted  once  to  the  combined  data  from  all 
individuals;  then  Information  on  Individual  variation  Is  represented  by  the 
distribution  of  observed  values  at  given  times  about  the  values  predicted  by 
the  model.  A  statistic  summarizing  this  variation,  and  potentially  useful  In 
describing  individual  variation,  is  the  mean  square  error,  which  is 
generally  reported  by  least-squares  optimization  routines. 

Uncertainties  in  parameters  fitted  simultaneously  by  an  optimization 
program  are  not  independent,  since  for  a  given  vector  of  fitted  parameters  a 
change  in  the  value  of  one  parameter  can  be  compensated  for  by  changes  in  the 
other  fitted  parameters.  For  example,  since  low-dose,  first-order  metabolic 
rate  is  approximated  by  vmax/km,  it  may  be  the  case  that  among  combinations  of 
vmax  and  km  providing  similar  fits  to  a  given  data  set,  a  positive  correlation 
exists  between  the  values  of  vmax  and  km.  Such  a  relationship  can  be  modeled 
approximately  by  assuming  a  bivariate  lognormal  distribution  for  vmax  and  km 
with  an  appropriate  correlation  parameter.  If  the  fitting  routine  employed  in 
an  optimization  provides  an  estimated  covariance  matrix  for  the  parameter 
estimates,  this  matrix  is  potentially  useful  in  estimating  a  value  of  the 
correlation  parameter.  At  this  time  we  have  no  basis  for  estimating  the 
correlation  between  vmax  and  km;  they  have  tentatively  been  treated  as  if  they 
were  independent. 

Estimation  of  Uncertainty  for  PERC  Circulatory  and  Ventilatory  Flow 
Rates .  This  category  of  par£imeters  included  ventilation  rate  (qp)  and  cardiac 
output  (qc)  (both  in  liters/hour) ,  and  percentages  of  cardiac  output  directed 
to  specific  compartments:  qlc  to  liver,  qrc  to  "rapidly  perfused,"  qsc  to 
"slowly  perfused,"  and  qfc  to  fat.  (In  some  cases  it  was  more  convenient  to 
treat  absolute  compartment  flows,  in  llters/Tir:  ql  for  liver,  qr  for  rapidly 
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perfused,  etc.)  Because  for  physiological  reasons  these  parameters  are 
correlated  (and  also  probably  correlated  with  compartment  volumes),  It  was 
convenient  to  treat  them  together. 

A  complete  accounting  of  uncertainty  might  Involve  detailed  assumptions 
regarding  physiological  mechanisms  regulating  the  relationships  among  the 
parameters.  The  description  here  Is  sufficiently  complete  to  demonstrate  the 
sort  of  physiological  considerations  that  are  relevant,  along  with 
possibilities  for  quantifying  the  way  that  these  considerations  may  determine 
uncertainty . 

The  formal  description  of  uncertainty  involved  the  three  following 
assumptions.  (1)  Uncertainty  with  respect  to  compartment  flows  (ql,  qr,  etc.) 
was  related  to  uncertainty  with  respect  to  corresponding  compartment  volumes. 
Volume  -  specif ic  perfusion  rates  (e.g.,  liters  blood/hr/liters  tissue)  were 
assumed  to  be  constant  so  that  the  rate  of  delivery  of  blood  to  a  compartment 
varies  as  the  compartment  volume  varied.  (2)  Additional  uncertainty  was 
incorporated  for  the  flows  in  humans.  This  added  uncertainty  related 
primarily  to  variation  in  the  levels  of  physical  exertion  in  daily  activities 
and  primarily  affected  qc,  qp ,  qs ,  and  qf.  (3)  Ventilation  rate  (qp)  was 
related  linearly  to  cardiac  flow  rate  (qc). 

Support  for  assumption  1  came  from  consideration  of  the  relationship 
between  individual  variation  in  compartment  volumes  and  individual  variation 
in  corresponding  blood  flow  rates.  Increase  in  a  compartment  volume  should  be 
accompanied  by  an  increase  in  blood  supply,  depending  on  the  relationship 
between  compartment  volume  and  oxygen  demand  per  unit  volume.  We  had  no  data 
relevant  to  the  latter  relationship,  and  so  we  assumed  constant  weight- 
specific  perfusion  rates  for  each  compartment.  Some  uncertainty  is  associated 
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with  the  values  of  the  weight  -  specif Ic  perfusion  rates.  We  did  not  quantify 
this  uncertainty;  we  assumed  that  perfusion  rates  were  constant  and  were  given 
by  the  ratios  of  preferred  compartment  flows  (in  liters  blood/hr)  to 
corresponding  preferred  compartment  volumes  (in  liters  tissue). 

The  assumption  of  constant  perfusion  rates  accounted  for  some  individual 
variation  in  compartment  flow  rates.  That  is,  as  compartment  volumes  varied, 
according  to  the  scheme  for  incorporating  uncertainty  in  the  volume  estimates 
(see  the  discussion  of  the  Dirlchlet  distribution  modeling  of  volume 
uncertainties  given  above),  the  flows  varied,  but  the  relationship  between 
flow  rate  and  compartment  volume  was  maintained. 

Additional  uncertainty  for  humans  with  given  compartment  volumes  may 
relate  primarily  to  the  level  of  physical  exertion  required  in  their  daily 
activities.  Increasing  activity  levels  result  in  higher  values  of  qc  and 
therefore  increasing  values  for  one  or  more  compartment  flows.  There  is 
evidence  that  the  increased  cardiac  output  is  directed  almost  entirely  to  the 
muscle.  Folkow  and  Neil  (1971)  reported  that  a  change  in  activity  level  from 
resting  to  heavy  exercise  resulted  in  approximately  equivalent  increases  in  qc 
and  flow  to  muscle  (qs) ,  for  average  subjects  and  top  athletes.  Data 
summarized  by  Hattis  et  al .  (1986)  were  roughly  consistent  with  this 
conclusion  (Table  VI-1-27):  the  discrepancies  may  be  due  to  the  fact  that 
different  entries  in  that  table  were  taken  from  different  primary  sources. 

Hattis  et  al .  (1986)  assumed  that  perfusion  of  fat  and  muscle  was 
related  to  activity  level,  and  that  perfusion  of  other  tissues  was  independent 
of  activity  level.  (The  conclusion  that  fat  perfusion  is  related  to  activity 
level  is  based  on  Astrand  (1983).)  Since  fat  is  typically  a  smaller 
proportion  of  body  weight  than  muscle,  this  does  not  contradict  the  conclusion 
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that  activity-related  increase  in  cardiac  output  is  allocated  primarily  to 
muscle.  We  assumed  simply  that  for  an  activity- related  increase  in  qc  of  dqc 
liters/hour,  qs  increased  by  dqs  -  dqc*vs/(vs+vf )  and  qf  increased  by  dqc-dqs, 
where  vs  was  the  volume  assumed  for  the  slowly  perfused  (muscle)  compartment 
and  vf  was  the  volume  assumed  for  fat.  It  was  convenient  to  assume  that  dqc 
was  distributed  lognormally  with  log-scale  mean  equal  to  zero  and  with  an 
uncertainty  factor  to  be  estimated. 

Regarding  the  third  assumption,  some  relationship  between  qp  and  qc  is 
to  be  expected  on  physiological  grounds.  Both  parameters  are  related  to  the 
oxygen  demand  of  tissues,  and  both  are  reported  to  follow  an  allometrlc 
relationship  with  body  weight  across  species,  with  exponent  about  0.8  (Stahl, 
1967).  Reitz  and  Nolan  (1986)  assumed  that  the  two  parameters  were  equal  in 
all  species;  data  from  various  sources  summarized  by  Hattis  et  al.  (1986)  (cf. 
Table  VI- 1-27)  suggested  that  such  an  equality  may  hold  approximately  for 
individuals  at  rest  but  that  with  Increasing  activity  levels,  qp  Increases 
more  rapidly  than  does  qc,  so  that  the  parameters  differ  by  a  factor  of  1.6 
with  light  exercise  and  by  a  factor  of  2.4  with  heavy  exercise.  Regarding  the 
linearity  of  the  relationship,  the  data  in  Table  VI -1-27  suggested  that  the 
true  relationship  over  activity  levels  may  be  convex:  as  activity  level 
changed  from  resting  to  "light  exercise,"  qp  increased  by  540  liters/hour 
while  qc  Increased  by  180  liters/hour,  a  ratio  of  Increments  of  3;  going  from 
rest  to  heavy  exercise,  the  ratio  of  increments  had  a  value  of  4.  However, 
the  range  of  activity  levels  represented  in  Table  VI-1-27  appears  to  be  much 
wider  than  the  range  of  activity  levels  characteristic  of  usual  daily 
activities,  which  might  be  generally  between  sitting  and  light  exercise.  In 
this  range  the  assumption  of  a  linear  relationship  may  provide  a  good 
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approximation.  We  approximated  the  relationship  by  a  line  with  slope  3 
connecting  the  "sitting"  values  (qp  -  420;  qc  -  440)  to  the  "light  exercise" 
values  (qp  -  960;  qc  -  582).  .\n  expression  of  the  assumption  that  was 

convenient  for  application  in  different  species  is 

(19)  qp  -  qpo  -  3  *  (qc  -  qco) , 


where  qpo  and  qco  are  the  "resting  values"  of  qp  and  qc ,  respectively,  for  a 
given  species.  For  rodents  we  tentatively  equated  resting  values  with 
preferred  values  in  Reitz  and  Nolan  (1986);  for  humans  we  used  the  "sitting" 
values  from  Table  VI-1-27. 

A  summary  of  the  uncertainty  considerations  associated  with  the  flow 
parameters  is  provided  here.  This  is  expressed  in  terms  of  the  algorithm  that 
was  implemented  to  generate  flow  parameter  values  (as  well  as  other  parameter 
values,  especially  the  compartment  volumes)  when  randomly  selected  parameter 
values  were  needed  as  input  for  dose  surrogate  uncertainty  estimation.  Note 
that  the  volumes  were  generated  by  the  appropriate  Dirichlet  distribution,  as 
discussed  above,  and  the  algorithm  then  maintained  the  relationship  between 
volume  and  flow  rate  discussed  in  the  preceding  paragraphs.  The  steps  taken 
were : 

1.  Generate  compartment  volumes  from  the  Dirichlet  distribution 
describing  their  uncertainty.  Denote  their  values  by  vi  (1  -  l,r,s,f). 

2.  Calculate  a  compartment  flow  rate,  qi,  for  each  compartment 
(1  -*  l,r,s,f)  according  to  the  equation 
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(20)  qi  -  Vi  *  (qip/vlp)  . 


where  qlp  and  vip  are  the  preferred  values  for  compartment  i  flow  rate  and 
volume,  respectively.  Calculate  qc  —  ql  +  qr  +  qs  +  qf. 

3.  For  humans  only.  Randomly  select  the  value  of  dqc  from  Its  lognormal 
uncertainty  distribution  (with  mean  zero  and  estimated  uncertainty  factor)  and 
increment  cardiac  output  and  flow  rates  to  the  fat  and  slowly  perfused 
tissues ; 


(21a) 

qc  -  qc  +  dqc ; 

(21b) 

qs  -  qs  +  dqc*(vs/(vs+vf ) ) ; 

(21c) 

qf  -  qf  +  dqc*(vf/(vs+vf ) ) . 

4.  Recalculate  percentage  of  cardiac  output  flowing  to  the  compartments: 


(22)  qic  -  qi/qc,  i  -  l,r,s,f. 


5.  Compute  ventilation  rate  according  to  formula  (19). 

In  the  Implementation  of  this  procedure,  preferred  values  for  flows  and 
volumes  were  equated  to  the  values  suggested  by  Reitz  and  Nolan  (1986)  with 
the  exception  that,  for  humans,  the  "shoeworker"  values  from  Table  VI-1-27 
were  used.  A  value  for  the  uncertainty  factor  associated  with  dqc  was  also 
required.  It  was  assumed  that  usual  daily  activities  vary  between  resting  and 
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light  exercise,  so  the  qc  values  associated  with  those  levels  of  activity  were 
taken  as  upper  and  lower  bounds  for  a  range  of  acceptable  values.  Following 
the  suggestions  given  above  for  estimation  of  uncertainty  factors  (cf.  Table 
VI-1-17),  the  uncertainty  factor  was  equated  to  statistic  SI,  which  in  this 
case  took  the  value  (582/400)°  ^  -  1.2,  where  582  and  400  are  the  appropriate 
qc  values  from  Table  VI-1-27. 

2.  Propagation  of  Parameter  Uncertainties:  Dose  Surrogate  and  Risk 
Uncertainty 

Having  estimated  uncertainty  distributions  appropriate  for  all  the  PBPK 
model  parameters  (for  mice,  rats,  and  humans)  one  can  determine  the  effects  of 
the  parameter  uncertainties  on  estimates  of  dose  surrogates  and  risks.  Monte 
Carlo  simulation  provides  the  basis  for  estimating  dose  surrogate  and  risk 
uncertainties . 

Appendix  VI-l-A  is  a  reprint  of  a  published  document  describing  an 
example  of  this  approach  to  risk  uncertainty  estimation.  The  example  extends 
the  work  presented  above  regarding  estimation  of  PERC  model  parameter 
uncertainties.  Bioassay  results  for  female  mice  exposed  to  PERC  are  used  for 
illustration. 
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Table  VI- 1-1 


Parameter  Sets  Used  In  Sensitivity  Analyses 


Parameter 

Mouse 

Rat 

Human 

bw 

0.025 

0.25 

70.0 

tiv 

0.003 

0.003 

0.003 

tgav 

vmaxc 

3.6 

1.27 

0.346 

km 

2.22 

4.56 

4.56 

pb 

24.43 

18.8 

10.3 

pi 

4.05 

3.74 

5.88 

Pf 

96.6 

87.1 

119.0 

ps 

5.2 

1.06 

3.10 

pr 

4.61 

3.74 

5.88 

vie 

0.04 

0.04 

0.0314 

vfc 

0.05 

0.07 

0.231 

vse 

0.78 

0.75 

0.62 

vre 

0.05 

0.05 

0.0371 

qpc 

28.0 

15.0 

15.0 

qcc 

28.0 

15.0 

15.0 

qlc 

0.25 

0.25 

0.24 

qfc 

0.05 

0.05 

0.05 

qsc 

0.19 

0.19 

0.19 

qrc 

0.51 

0.51 

0.52 
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Table  VI-1-2 


Doses  Employed  In  Sensitivity  Analyses 


Dose  Pattern 

Low 

Medium 

Hi  eh 

Gavage 

Mouse 

0.30 

30 

2000 

(doses  in  mg/kg) 

Rat 

0.15 

15 

1400 

Human 

0.05 

6.0 

600 

Intravenous  injection 

Mouse 

0.01 

2.0 

200 

(doses  in  mg/kg) 

Rat 

0.03 

3.0 

350 

Human 

0.04 

4.0 

450 

Inhalation,  8  hrs/day 

Mouse 

0.30 

30 

1600 

(doses  in  ppm) 

Rat 

0.60 

60 

5000 

Human 

1.0 

115 

11500 

Inhalation,  continuous 

Mouse 

0.30 

(24) 

24  (24) 

1400  (48) 

(doses  in  ppm*) 

Rat 

0.50 

(48) 

44  (48) 

3800  (48) 

Human 

0.90 

(24) 

90  (24) 

9000  (48) 

•In  parentheses  are  the  times  (in  hours)  that  mark  the  beginning  of  the 
24-hour  period  used  to  calculated  daily  dose  surrogates.  The  times  are 
those  at  which  the  system  is  close  to  steady  state. 
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TabU  VI-l-3 
pate  Hodal  Smaltlvlty 
Houa*  ParsMtar  Sat:  Oavasa  Adalnlatration 
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Sanaltlvlty  la  aMpraaaad  oa  pareant  chonga  in  ttia  doaa  aurrogota  induced  by  o  1<  chonge  in  the  Indlcoted  poronatar. 
Tha  oatarlaka  algnlfy  on  obaoluta  mognltuda  laaa  tt>on  I.OOE-5. 
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Sansltlvlty  l(  •Kpraaawf  os  psrcsnt  chongs  in  tho  doss  surrogots  Inducsd  by  a  l|f  cbongs  In  tbs  Indicated  poratnstsr. 
Ths  ostsrlsks  signify  on  absolute  mognltude  loss  than  1.00E-b. 
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PERC  Hod*l  Sansltlvlty 

Moua*  Parainatar  Sat:  8  Hour  Inhalation  Athalnlatratlon 
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Senaltltflty  la  eaproased  oa  psrcent  change*  In  the  dose  surrogote  Induced  by  o  l|(  chonge  in  the  Indlcoted  porometer. 
The  oaterlalia  algnlfy  on  obaolute  mognltude  lesa  thon  1.00C-S. 
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PEXC  Hodal  S«naltlvlty 

Houaa  Paramatar  Sat;  Contlnuoua  Inhalation  Adalnlatratlon 
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PfXC  Modal  Sanaltlvlty 
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Sanaltlalty  la  aapraaaad  ea  pareant  ehonga  in  tha  doaa  aurrogota  Inducad  by  o  ehongo  In  tha  indleotad  poromatar. 
Tha  aatarlaha  algnlfy  an  obaaluta  laognltuda  laaa  than  l.OOC-0. 


Table  VI-l-lS 


Relative  Ranking  of  Paraiueter  Sensitivity  of 
Virtual  Concentration  of  Metabolite;  Mouse  Parameter  Set 


vie-* 

vlc- 

vie  - 

vlc  • 

vlc- 

vmaxc+ 

qpc- 

qpc- 

pb+ 

qpc+ 

vmaxc+ 

vlc- 

pb+ 

pb+ 

qpc- 

pb+ 

pb+ 

bw- 

lan- 

viiiaxc  + 

vTiiaxc+ 

vmaxc+ 

ko- 

\niiaxc+ 

km- 

kn- 

km- 

bw- 

qcc  - 

pf+ 

vfc+ 

bv- 

qpc+ 

ps+ 

VSC  + 

qcc+ 

qlc+ 

qcc+ 

qcc+ 

qcc+ 

qlc+ 

qsc+ 

qlc+ 

qlc+ 

qlc+ 

pl+ 

pr+ 

pr+ 

vrc+ 

vrc+ 

qfc+ 

qfc+ 

VSC  + 

pl+ 

ka- 

ps+ 

pb+ 

pf+ 

vfc+ 

pl+ 

qrc- 

qsc- 

pf. 

qpc+ 

qsc+ 

bw- 

bw- 

vfc- 

vfc* 

qfc+ 

vrc+ 


qrc- 

pf- 

ps- 

qcc+ 

vfc- 

vse- 

qlc+ 

qfc+ 

qfc+ 
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Table  VI-1-15  (continued) 


Relative  Ranking  of  Parameter  Sensitivity  of 
Virtual  Concentration  of  Metabolite;  Mouse  Parameter  Set 


Range  of 
Percent  Changes 


_ IV  Dosing _  Continuous  Inhalation  Exposure 

L _ M _ H _ L _ M _ H 


ps- 

qsc+ 

pf- 

vsc  - 

pr- 

vfc- 

qsc+ 

vrc- 

ps- 

vrc+ 

pl- 

vsc  - 

pl+ 

qrc  + 

qfc+ 

pr+ 

qsc+ 

qrc  + 

vrc- 

qrc* 

pl+ 

pr+ 

^The  plus  or  minus  sign  following  the  parameter  indicates  the  direction  of 
change  of  CM  given  1%  increase  in  the  parameter  value. 
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Table  VI. 1-16 


Range  of 
Percent  Changes 


►  lOO 


10-1-  -  lOO 


10-2  .  10-1 


10-3  .  10-2 


10-4  .  10-3 


10-5  .  10-4 


L 


aclve 

Ranking  of 

Parameter 

Sensitivity  of 

Concentration  of 

He  ta^'ol  i  t 

e;  Rat  Pai 

aeter  Set 

IV  Dosing 

Continuous 

Inhalation  Exoosure 

L 

M 

H 

L 

H 

H 

vniaxc+* 

vie  - 

vie- 

vlc- 

vlc- 

vlc- 

vlc- 

vniaxc+  pb+ 

\Tiiaxc+ 

pb+ 

vniaxc+ 

bw- 

km- 

qpc- 

qpc- 

viDaxc  + 

pb+ 

pb+ 

viiiaxc+- 

pb+ 

ko- 

kn- 

qpc- 

ka- 

km- 

bw- 

bw- 

qcc+ 

qcc+ 

qpc+ 

qpc+ 

vfc  - 

qfc  + 

pf- 

bw- 

vfc  - 

qfc-*- 

VSC  + 

qcc+ 

qcc+ 

pf- 

qlc+ 

ps+ 

qlc+ 

qlc+ 

qcc+ 

bw- 

vfc+ 

qlc+ 

pl+ 

pf+ 

qfc+ 

VSC  + 

qlc+ 

bw- 

ps  + 

qsc+ 

pr+ 

vrc+ 

ps- 

qrc  + 

vrc+ 

vfc- 

pf- 

kn- 

vse  - 

pr+ 

pf- 

vfc- 

pb+ 

pl+ 

qfc+ 

qfc+ 

qsc 

pr- 

vse- 

vse- 

vfc- 

VTC- 

ps- 

ps- 

pf- 

qsc+ 

qpc  + 

pl- 

qrc+ 

qrc- 

VTC- 

pr- 

qcc+ 

pr- 

vre- 

qfc+ 

pl- 

pl- 
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Table  VI-1-16  (continued) 

Relative  Ranking  of  Parameter  Sensitivity  of 
Virtual  Concentration  of  Metabolite;  Rat  Parameter  Set 


Range  of  _ IV  Dosing  _ _  Continuous  Inhalation  Exposure 

Percent  Chanees  _ L _ M  H _ L _ M _ __H _ _ 

qsc+  qsc+  vse- 

qrc  qrc+  ps- 

qlc+ 
vre- 
pr- 
pl- 
qsc+ 
qrc 


^The  plus  or  minus  sign  following  the  parameter  Indicates  the  direction  of 
change  of  CM  given  1%  increase  in  the  parameter  value. 
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Table  VI-1-17 


Statistics  Relevant  to  Uncertainty  Estimation 


GM 


Geometric  mean:  antilog  of  mean  log. 


SSDL 


CV 


SI  -  (UB/LB)0-5 


S2  -  exp(2*SDLE) 


Sample  standard  deviation  of  natural  logs. 

The  maximum  likelihood  and  preferred  estimator 
for  SDL,  the  log-scale  standard  deviation. 

Sample  coefficient  of  variation 
-  standard  deviation/mean.  Used  here  often 
to  approximate  the  log-scale  standard 
deviation. 

Use  as  an  estimate  of  a  UF  based  on  a 
reference  Interval  LB  to  UB  (usually 
representing  expert  opinion.) 

for  SOLE  an  estimate  of  the  log- scale  standard 
deviation,  (either  SSDL  or  CV.)  Use  as  a 
direct  estimate  of  a  UF  for  an  irreducible 
uncertainty  source  based  on  data 
assumed  ild*. 


S3  -  exp(2*SDLE/SQRT(#M))  When  the  prefered  value  is  the  GM  of  iid 

measurements,  use  S3  as  an  estimator  for 
the  UF. 


®  iid  -  independent,  identically  distributed, 
from  the  same  probability  distribution.) 


(Generated  independently 
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Table  VI-1-18 


Data  Corr.plled  by  Fi  s  e  r  ova  -  Be  r  ge  rove  (1983)  Regarding  Study- to-Study 
Variation  in  Estimates  of  Blood/Gas  Partition  Coefficients  at  37°C 


Chemical 


Species  Estimate  (Reference®) 


Acetone 

roan 

Acetone 

dog 

Benzene 

roan 

Carbon  Tetrachloride 

man 

Chloroform 

roan 

Cyclopropane 

roan 

1 , 1 -Dichloroethane 

man 

1 , 2 -Dichloroethane 

man 

Diethylether 

man 

Diethylether 

dog 

Methylene  Chloride 

roan 

Ethane 

roan 

Ethane 

dog 

Fluroxene 

roan 

Halothane 

man 

Halothane 

dog 

1,1,2 , 2-Te trachloro- 

roan 

ethane 

PERC 

man 

Toluene 

roan 

S  tyrene 

roan 

Trichloroethylene 

man 

2^.5(l),  3^1(2),  313(3).  302(4.) 
363(2),  376(3).  291(6) 

7.8(1).  6.5(5),  9.0(6). 

7.8(7).  7.7(8) 

2.4(2).  0.6(9),  5.8(6) 

10.3(10).  9.0(9).  8.2(5).  11.0(6) 
0.58(2).  0.35(3),  0.43(11).  0.51(4) 
4.7(10).  4.5(9).  6.0(6) 

19.5(10),  20(9).  19.7(6) 

12.5(2),  11.6(3).  12.3(4), 

12.8(6),  11.9(11) 

12.7(2),  10.7(4).  11.8(3) 

9.7(10).  7.0(9),  9.4(6) 

0.1(2),  0.05(3),  0.08(4) 

0.13(2).  0.04(3).  0.10(4) 

1.3(12),  1.5(2),  1.4(11) 

2  6(12).  2.5(2),  2.6(13), 

2.8(3),  2.5(11).  2.5(14) 

1.7(3),  3.5(38),  3.8(2) 

121(10),  73(9),  141(6) 

13.1(10),  9.1(9),  18.9(6) 

15.6(1),  15.6(5),  14.7(8),  16.3(6) 
32(30),  52(1),  64(6) 

9(5),  9.5(9).  9.9(6),  9.9(15) 


®(1)  Sato  and  Nakajima  (1979);  (2)  Uagrier  et  «J..  (1974);  (3)  Young  and 
Vagner  (1979);  (4)  Dueck  et  al..  (1978);  (5)  Sherwood  (1976);  (6)  Bocek 
(personal  communication  to  Flserova-Bergerova ,  ibid.);  (7)  Telslnger 
and  Skramovsky  (1947);  (8)  Sato  eC  ml..  (1972);  (9)  Morgan  et  si.. 
(1970);  (10)  Sato  and  Makajlma  (1979);  (11)  Gibbs  et  al . .  (1975);  (12) 
Ellis  and  Stoeltlng  (1975);  (13)  Saraiva  et  al . .  (1977);  (14)  Eger 
et  si.  (1962);  (15)  Fink  and  Morlkawa  (1970). 
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Table  VI-1-19 


Analysis  of  Inter-Study  Variation  In  Estimates  of  Blood/Alr 
Partition  Coefficients:  Summary  of  Data  in  Table  VI-1-18 
(From  Data  Compiled  by  Flserova-Bergerova  (1983)) 


Chemical  (Species) 

GM^ 

SSDL^ 

S2C 

^Studies 

Ethane 

(man) 

0.074 

0.35 

2.0 

3 

Ethane 

(dog) 

0.080 

0.62 

3.4 

3 

Cyclopropane 

(man) 

0.46 

0.22 

1.5 

4 

Fluroxene 

(man) 

1.4 

0.072 

1.2 

3 

Carbon  Tetrachloride 

(man) 

2.0 

1.1 

9.8 

3 

Halothane 

(man) 

2.6 

0.044 

1.1 

6 

Halothane 

(dog) 

2.8 

0.44 

2.4 

3 

1,1-Dlchloroe thane 

(man) 

5.0 

0.16 

1.4 

3 

Benzene 

(man) 

7.7 

0.12 

1.3 

5 

Methylene  Chloride 

(man) 

8.6 

0.18 

1.4 

3 

Chloroform 

(man) 

9.6 

0.13 

1.3 

4 

Trichloroethylene 

(man) 

9.6 

0.045 

1.1 

4 

Diethylether 

(dog) 

12 

0.086 

1.2 

3 

Diethylethe  r 

(man) 

12 

0.039 

1.1 

5 

*  PERC 

(man)^ 

12 

0.32 

1.9 

4 

Toluene 

(man) 

16 

0.042 

1.1 

4 

1 , 2 • D1 chi o roe  thane 

(man) 

20 

0.013 

1.0 

3 

Styrene 

(man) 

47 

0.36 

2.0 

3 

1.1,2,2-Tetra- 
chloroe  thane 

(man) 

108 

0.34 

2.0 

3 

Acetone 

(roan) 

298 

0.14 

1.3 

4 

Acetone 

(dog) 

341 

0.14 

1.3 

3 

^Geometric  mean  of  available  estimates.  The  table  is  sorted  on  this 
column,  to  demonstrate  any  relationship  betveen  CM  and  SSDL. 

^Standard  deviation  (In  In  scale)  of  available  estimates. 

CS2  -  exp(2*SSDL). 

^In  addition  to  the  values  reported  in  Table  IV-1-18,  a  value  of  10.3  has 
been  Included  based  on  communication  of  M.  Andersen  to  Reitz  and  Nolan 
(1986). 
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Tabl*  VI-1-20 


Estimates  of  the  Log-Scale  Standard  Deviation  (SOLE)  Representing 
Variation  Among  Individual  Humans  and  Animals  in  Tissue/Alr  Partition 
Coefficients,  for  Trichloroethylene,  Haloethane,  Cyclopropane  and  NO2 


Trichloroethylene*  Halothane* 

Halothane^ 

5  Humans 

5  humans 

18-38  horses 

Mean  SDLE 

Kean  SDLE 

Mean  SDLE 

Blood 

8 

0.3 

2 

0.2 

2 

0.2 

Brain 

21 

0.2 

6 

0.5 

5 

0.1 

Kidney 

15 

0.3 

6 

0.3 

3 

0.4 

Liver 

29 

0.4 

9 

0.5 

9 

1.0 

Lung 

14 

1.0 

4 

0.7 

3 

0.2 

Muscle 

19 

0.3 

8 

0.2 

4 

1.0 

Fat 

569 

0.03 

222 

0.  1 

Cyclopropane^ 

NO2 

c 

5-9  Rabbits 

5-6  Rabbits 

Max 

SDLEd 

Mean 

SDLE 

Mean 

SDLE 

Blood 

0.7 

0.04 

0.5 

0.01 

0.3 

Brain 

r  .8 

0.1 

0.4 

0.02 

0.5 

Kidney 

0.9 

0.1 

0.4 

0.04 

0.4 

Liver 

0.7 

0.1 

0.4 

0.03 

1.0 

Lung 

1.0 

Muscle 

0.51 

0.1 

0.4 

0.07 

1.0 

Fat 

0.1 

^SDLE  CV,  computed  from  standard  deviations  and  means  reported  by 
Fiserova-Bergerova  et  si.  (1984),  each  based  on  five  subject-specific 
estimates  (the  sane  five  subjects  for  each  chemical).  Tissues  were 
from  autopsied  Individuals. 

^SDLE  from  Vcbb  and  Vcaver  (1981),  each  the  square-root  of  a  variance 
component  computed  in  the  log  scale. 

CSDL£  based  on  variance  components  reported  by  Mapleson  et  si.  (1970). 
Since  computations  were  carried  out  in  the  log  scale,  square  roots  of 
reported  variance  components  are  estimates  of  the  log- scale  standard 
deviations.  For  NO2 ,  ve  use  the  reported  betveen-anlmals  component; 
for  cyclopropane  a  separate  betveen-animals  component  was  not 
reported,  and  so  the  estimate  used  incorporates  relatively  more 
measurement  error.  For  cyclopropane,  our  CV  is  the  square  root  of  the 
average  of  the  'la  vivo*  and  "tn  vitro"  variance  component  estimates. 
Multiple  strains  of  rabbits  were  used  In  the  experiment. 

^SDLE  for  a  tissue  is  the  maximum  of  the  SDLE  tabulated,  over  species 
and  studies. 
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Table  VI-1-21 


Specles-to-Specles  Variation  In  Measurements  of  Tlssue/Air 
Partition  Coefficients  (Flserova-Bergerova,  1983). 


1 SOFLURANE^ 


Tissue 

1 

1 

_ 1 _ 

Man 

Monkey 

Dog 

Rat 

1 

1 

1 

GMt> 

Brain 

1 

1 

1.9 

2.0 

2.5 

2.2 

1 

1 

2.1 

Heart 

1 

1.8 

1.6 

2.7 

2.8 

1 

2.2 

Kidney 

1 

2.1 

1.9 

3.6 

2.5 

1 

2.4 

Liver 

1 

4.1 

4.6 

4.2 

3.2 

1 

4.0 

Lung 

1 

1.6 

1.5 

2.6 

2.4 

1 

2.0 

Muscle 

1 

2.4 

1.4 

3.4 

1.6 

1 

2.1 

Fat 

1 

_ 1 _ 

69 

66 

75 

63 

1 

_i_ 

68 

METHYLENE  CHLORIDgc 


Tissue 

Man 

Monkey 

Dog 

Rat 

i  GM= 

Brain 

7.1 

6.8 

7.9 

7.4 

1  7.3 

Heart 

7.1 

6.3 

7.4 

9.2 

1  7.4 

Kidney 

5.8 

6.3 

9.5 

7.4 

1  7.1 

Liver 

7.1 

11 

11 

8.9 

1  9.4 

Lung 

5.8 

5.5 

7.6 

7.9 

1  6.6 

Muscle 

1  ^.7 

3.7 

7.9 

4.7 

1  5.0 

Fat 

1  84 

1 _ 

86 

97 

91 

1  89 

®Data  from  Table  3  of  Fiscrova-Bergerova,  1983. 
^Geometric  mean  of  four  specles-speclflc  estimates. 
^Estimated  from  Figure  5  of  Flserova-Bergerova ,  1983. 
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Table  VI-1-22 


Analysis  of  Spec ies - to- Spec les  Variation  In  Tlssue/Gas  Partition 
Coefficients,  for  Isofluorane  and  Methylene  Chloride 
(Data  Compiled  by  Flserova-Bergerova,  1983; 

Species  are  Man,  Monkey,  Dog.  and  Rat.) 


ISOFLURANE  METHYLENE  CHLORIDE 

Human/Rat  Human/Rat 


Tissue 

1 

J _ 

CM 

Error^ 

S2 

1 

_  1 

CM 

Error* 

S2 

1 

Brain 

1 

1 

1 

2.1 

1.2 

1.3 

1 

1 

1 

7.3 

1.0 

1.1 

1 

1 

Heart 

I 

1 

1 

2.2 

1.6 

1.7 

1 

1 

1 

7.4 

1.3 

1.4 

1 

1 

Kidney 

1 

1 

1 

2.4 

1.2 

1.8 

1 

1 

1 

7.1 

1.3 

1.6 

! 

1 

Liver 

1 

1 

1 

4.0 

1.3 

1.3 

1 

1 

9.4 

1.3 

1.5 

1 

1 

Lung 

1 

1 

1 

2.0 

1.5 

1.7 

1 

1 

6.6 

1.4 

1.4 

1 

1 

Muscle 

1 

1 

1 

2.1 

1.5 

2.3 

1 

1 

5.0 

1.0 

2.0 

1 

1 

1 

Fat 

1 

1 

J _ 

68 

1.1 

1.2 

1 

1 

1 

89 

1.1 

1.1 

1 

1 

^Let  R  be 

the 

ratio  of 

the  human 

value 

to  the 

rat 

value;  the 

value 

reported  is  the  larger  of  R  and  1/R,  (This  corresponds  to  taking  the 
absolute  difference  in  the  log  scale,  then  transforming  back  to  the 
original  scale.) 
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Table  Vl-1.23 


SujTurary  of  PERC  Partition  Coefficient  Uncertainty  Estimation 


Source  M: 

1.  Estimate  from  data  SOLE  s  for  directly  measured  tissue/alr  PC's; 

a.  SMl(l.b)«  1  -  r.m.h. 

b.  SMKr.l) 

c.  SMl(r.m) 

d.  SMl(r.f) 

2.  Estimate  from  data  SDLE's  corresponding  to  variations  across  species  of 
tissue/air  PC’s;  SM2(J)b,  J  -  b.l.r.n.f. 

3.  Estimate  SDLE  corresponding  to  difference  between  log  of  liver/air  PC 
and  log  of  rapid/alr  PC;  SM3(l,r)^. 

h.  Compute  SDLE's  for  human  or  mouse  tlssue/air  PC's  (tissues  other  than 
blood  or  rapidly  perfused  organs)  when  these  are  based  on  measurements 
in  rats; 

SM^(i.j)  -  (SMl(r.j)2  +  2*SM2(j)2)0.5^  i  -  a.h;  j  -  l.m.f. 

5.  Compute  rapid/air  PC's  when  these  are  based  on  rat  llver/air 
measurements ; 

SM4(i.r)  -  [SMl(r.l)2  +  SH2(1)2  +  SM2(r)2  +  SM3(1 .r)2 jO.5^ 

1  -  r .m ,h . 

6.  Compute  SDLE's  and  UFj^'s  for  blood/air  and  tlssue/blood  PC's; 

a.  Blood/air,  for  species  i  (t  -  r.m.h)  - 

SDLECi.b}  -  SMl(l.b). 

b.  Tissue/blood,  for  rats  and  tissue  J  other  than  rapidly  perfused  (j  - 
l.m.f)  • 

SDLE(r.J)  -  [SMl(r.J)2  +  SMl(r.b)  )0- 5. 

c.  Tissue/blood,  for  species  I  other  than  rats  (i  -  m.h)  and  tissue  J 
(J  -  l,r,m,f)  or  for  rats  and  rapidly  perfused  tissue  (I  -  r, 

j  -  r)  - 

SDLE(i.J)  -  [SM4(i,j)2  +  SM1( i .b) 2  ]0 . 5 . 

In  each  case,  UFM(i,j)  -  exp(2*SDLE(l , J))  . 
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Table  VI-1-23  (continued) 

Su^wT.ary  of  PERC  Partition  Coefficient  Uncertainty  Estimation 


Source  1 : 

1.  Estimate  from  data  SDLE's  representing  Individual  variation  In 
tlssue/alr  PC's.  Denote  these  by 

SIl{i.J)‘*  1  -  r.m.h;  J  -  b.l.r.m.f. 

2.  Compute  SDLE's  and  UFi's  for  blood/air  and  tlssue/blood  PC's: 

a.  Blood/air,  for  species  i  (1  -  r.m.h)  - 

SDLE(l.b)  -  Sll(l.b). 

b.  Tissue/blood,  for  species  i  (1  -  r  m.h)  and  tissue  j  (j  "  l.r.m.f)  - 

SDLE(i.J)  -  (Sll(i.j)2  -t-  SI  1  (i  .b)  ^  ]  0 . 5  _ 

In  each  case,  Uri(i,j)  -  exp ( 2*SDLE( i  .  J ) )  . 

Total  Uncertainty: 

1.  Combine  source  M  and  source  I  uncertainty  factors  for  species  1  and 
tissue  j  (1  -  r.m.h;  J  -  b.l.r.m.f): 

UFjd.J)  -  exp[{ln(UFM(l,J))2  +  ln(UFi  ( 1 .  J  ) )  2  )0 . 5)  . 


“The  first  variable  In  parentheses  is  a  species  code  and  the  second  Is 
a  tissue  code.  This  notation  corresponds  to  that  in  Figure  Vl-1*2 , 
^The  variable  in  pa^'^ntheses  Is  a  tissue  code.  This  notation 
corresponds  to  that  in  Figure  VI-1*2. 

<^rhe  tvo  variables  in  parentheses  are  tissue  codes.  This  notation 
corresponds  to  that  in  Figure  VI -1-2 . 

<JThe  variables  In  parentheses  are  like  those  for  the  analogous  SMl 
standard  deviations,  i.e.,  a  species  code  followed  by  a  tissue  code. 
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Table  VI-1-24 


Source  1,  Source  M,  and  Total  Uncertainty  Factors 
for  PERC  Partition  Coefficient  Estimates 


PBPK  Parameter 

Preferred 

Parameter 

Value 

UFi 

UFm 

UFt 

RATS 

Blood/air 

18.8 

1.8 

2.0 

2.5 

Liver/blood 

3.74 

8.1 

2.7 

10 

Rapid/blood 

3.74 

3.2 

4.3 

6.5 

Huscle/blood 

1.06 

8.1 

2.7 

10 

Fat/blood 

87.1 

1.9 

2.7 

3.2 

MICE 

Blood/alr 

16.9 

1.8 

2.0 

2.5 

Liver/blood 

3.01 

8.1 

2.7 

10 

Rapid/blood 

3.01 

3.2 

4.3 

6.5 

Kuscle/blood 

2.59 

8.1 

2.7 

10 

Fat/blood 

48.3 

1.9 

2.7 

3.2 

HUKANS 

Blood/alr 

10.3 

1.8 

1.4 

2.0 

Live  r/blood 

5.88 

8.1 

3.1 

11 

Rapld/blood 

5.88 

3.2 

3.8 

5.8 

Muscle/blood 

3.10 

8.1 

4.1 

12 

Fat/blood 

119.1 

1.9 

3.1 

3.6 
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Table  VI-1-25 


Variation  In  Voiaxc  Estimates 


Source 

Es  t imates 

(mg/hrAgO-^^)  CM 

Reitz  & 
Nolan* 

SSDL^ 

Mice 

Hattls  eC  aJ.  "Schumann" 
Hattis  et  ai.  "Mitoma" 
Reitz  &  Solan  "Mouse  1" 
Reitz  6  Solan  "Mouse  2" 

1.53c  3.96 

5.34 

8 . 34 

3.60 

8.34 

0.721 

Rats 

Hattis  ec  al .  "Pegg" 

Hattis  et  al.  "Mitoma" 
Reitz  &  Nolan,  F344 

2,39  2.03 

2 . 78 

1.27 

1.27 

0.414 

Humans 

Hattis  ec  aJ.  "Ohitsu(+)" 
Hattis  et  al.  "Ik€d;(+)" 
Reitz  &  Nolan 

0.389  0.330 

0.266 

0.346<J 

0.346 

0.195 

Weighted  Average  SSDL* 

0.532 

^Preferred  values  from  Reitz  &  Nolan,  using  the 
mice . 

^SSDL  -  standard  deviation  of  In  estimates. 

^From  the  range  of  values  1.53-1.72  suggested  by 

".Mouse  1" 

Hattis  et 

estimates  for 

al . ,  the 

value  1.53  yields  the  largest  value  of  SSDL. 

value  of  0.256  Is  given  in  Table  1  of  our  draft  of  Reitz  and  Nolan. 
Ve  have  used  the  value  of  0.3^6  from  the  text  of  the  draft,  vhich  we 
believe  to  be  the  correct  value. 

'^Weighted  average  -  ( suni(df*SSDL2)/sua(df )  )0 ■  ^  where  sujnniation  is  over 
mice,  rats,  and  hurians,  and  df  is  the  number  of  measurements  available 
for  a  species,  minus  one  (e.g.  df  -  3  for  mice). 
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Table  Vl-1-26 


Variation  In  Km  Estimates 


Source 

Es  t imates 
(mg/llter) 

CM 

Reitz  & 
Nolan* 

SSDL^ 

MIcp 

Hattis  eC  si.  "Schumann* 
Hattis  et  si.  "Mitoma" 

Reitz  &  Nolan  "Mouse  2" 

0.268C 

5.35 

2.22 

1  .47 

4.56 

1 . 54 

Rats 

Hattis  et  al.  "Pegg" 

Hattis  et  al .  "Mitoma" 

Reitz  &  Nolan,  F344 

6.78 

17.8 

4.56 

8.19 

4.56 

0.700 

Humans 

Hattis  et  si.  "Ohitsu<+)" 
Hattis  et  si.  "Ikeda(+)" 

5.69 

0.  609 

1.86 

4.56 

1.58 

Weighted  Average  SSDL*^ 

1.28 

^Preferred  values  from  Reitz  &  Nolan,  equal  to  the  rat  estimate  for 
all  three  species. 

^SSDL  -  standard  deviation  of  In  estimates. 

^From  the  range  of  values  0.268-0.669  suggested  by  Hattis  et  a].,  the 
value  0.268  yields  the  largest  value  of  SSDL. 

'^Weighted  average  -  ( sum(df*SSDL2)/sun(df )  )0  •  5  where  summation  is  over 
mice,  rats,  and  humans,  and  df  is  the  number  of  measurements  available 
for  a  species,  minus  one  (e.g.  df  -  3  for  mice). 


Table  VI-1-27 


Hour  Volumes  (in  liters)  for  Ventilation  (qp)  ,  Cardiac  Output  (qc), 
and  Muscle  Blood  Supply  (qm)  ,  as  a  Function  of  Activity  Level  in  Humans 


(  f  rom 

Hattis  et 

al .  , 

1986) 

Activity  Level 

qp 

(incr . )* 

qc 

(  incr . ) * 

qm 

(incr. )* 

Sitting 

420 

400 

66 

" Shoeworker 

683 

(  263) 

489 

( 

89) 

Light  exercise 

960 

(  540) 

582 

( 

182) 

324 

(258) 

Heavy  exercise 

2040 

(1620) 

840 

( 

440) 

540 

(474) 

^Increment  from  sitting  value  -  value  at  specified  activity  level  minus 
corresponding  sitting  value. 

^qp  from  Brugnone  et  al.  (1980),  assumed  typical  for  occupational 
settings;  qc  estimated  by  linear  interpolation  betveen  adjacent  values 
of  qp  and  qc . 
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Evaluation  of  uncertainty  in  input  parameters  to 
pharmacokinetic  models  and  the  resulting 
uncertainty  in  output 
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(Acctplcd  10  July  1989) 

A>t  nords  L'nctriainiy  analysis,  PBPK  modeling.  Risk  assessmeni.  Teirachloroelhylene 


SUMMARY 

Pbysiolc'E’^'-'ly  based  pbarmacckinelic  (PBPK)  models  may  b<  used  lo  predict  the  concentrations  of 
parent  chem.-.l  or  metabolites  in  tissues,  resulting  from  specified  chemica-’  exposures  An  important  appb- 
cation  of  PBPK  modeling  is  in  assessment  ofi  carcinogenic  risks  to  humans,  based  on  arumal  data  The 
parameten  o{  a  PBPK  model  may  include  meubolic  p>atameters.  blood  air  and  tissue  blood  (sartition  eo- 
efficienls.  and  physiological  parameters,  such  as  organ  ueighu  and  blood  flow  rates  Uncertainty  in  esti¬ 
mates  or  these  parameters  results  in  uncertainty  regarding  tissue  concentnlioos  and  resulting  risks  Data 
are  re'iewed  relevant  to  the  quantification  ofi  these  uncertainties,  ror  a  PBPK  model-based  nsk  assessment 
Tor  teirachloroelhylene  Probability  distributions  art  developed  to  express  uncertainly  in  model  pararo- 
e;ers.  and  uncertainties  are  propagated  by  a  sequence  o(  operations  that  simulates  processes  recognized 
as  contributing  to  estimates  or  human  nsk  Distributions  of  PBPK  model  output  and  human  risk  estimates 
are  used  to  charactenze  uncertainty  resulting  from  uncertainty  in  model  parameten. 


INTRODUCTION 

Physiologically-bascd  pharmacokinetic  (PBPK)  models  can  be  used  in  the  animal- 
based  assessment  of  human  cancer  risks.  A  PBPK  model  is  assumed  for  rodents  and 
humans  (with  parameter  values  that  are  possibly  species-specific)  and  dose  surrogates 
art  calculated  on  the  basis  of  that  model.  A  dose  surrogate  is  a  particular  measure 
of  chemical  delivered  to  a  putative  target  tissue.  The  dose  surrogate  values  corre- 


Addrtu  for  (orrespondenct  B  Allen,  Clement  Associate},  Inc.,  1201  Caines  Street,  Rustoo,  LA  7I27Q, 
USA. 
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VI-l-8f 


spending  to  the  dose  levels  and  dose  regimen  follov^ed  in  a  chosen  carcinogenicity 
bioassay  and  the  observed  tumor  count  data  are  the  input  to  dose-response  model 
(e  g.,  the  multistage  or  Weibull  models  (1)).  Values  of  the  chosen  surrogate  dose  are 
estimated  for  humans,  corresponding  to  specified  external  exposures,  and  human 
risks  are  estimated  based  on  the  fitted  dose-response  curve. 

This  report  describes  the  development  of  methods  for  analyzing  the  effect  of  uncer¬ 
tainties  associated  with  PBPK  model  parameters.  The  methods  are  illustrated  in  the 
context  of  a  nsk  assessment  for  tetrachlorocthylene  that  employs  a  multicompart¬ 
ment  PBPK  model  and  a  carcinogenic  end-point  observed  in  female  mice  (2). 

PHYSIOLOGICALLY  Based  pharmacokinetic  model  for  tetrachloroeth- 

YUENE 


The  Environmental  Protection  Agency  (EPA)  (3]  and  Hattis  et  al.  [4]  have  devel¬ 
oped  two  slightly  different  PBPK  models  for  tetrachloroethylene.  Reitz  and  Nolan 
[5]  have  reviewed  the  parameter  values  adopted  ov  the  EPA  and  have  recommended 
some  revisions.  In  order  to  evaluate  uncertainties  within  the  framework  of  a  given 
PBPK  model,  analy  ses  presented  here  are  based  on  the  EPA  model  (Fig.  1).  Flow- 


Fig  I.  TelracKJorocili)leoe  PBPK  model.  Nouuoii:  €,  tod  e,  art  cooccnintioa  of  parent  in  in¬ 

haled  air,  exhaled  air.  arterial  blood,  and  venous  blood;  f.  r,  and  /,  are  Ihe  perfusioa  rale,  paical  coooca- 
tration,  and  tiuue>blood  paniiion  coeffident  for  comparUDCiii  c  and  and  are  constanu  detennia- 
ing  the  rate  of  meiatnlism  in  liver. 
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in 


limllcd  comparimenis  and  Michaclis-  Mcnlcn  metabolism  (occurring  in  the  liver)  are 
assumed. 

Notation  and  preferred  values  for  model  parameters  are  listed  in  Table  I.  In  gener- 


TABLE  I 

PARAMETERS  OF  THE  TETRACH LOROETHYLESE  PBPK  MODEL.  WITH  PREFERRED 
VALUES  AND  UNCERTAINTY  FACTORS(UF‘) 


Patamcler  Prcfcrrwl  valuti  (UFi) 


Mice  Humans 


Bod)  weight  (i»,  kg) 

0  028 

70 

Companmenl  proportions  (range  O-I) 

Liver  <v^) 

0  056  0  24) 

0  026 

(1.35) 

Rapioi)  perf  (v„) 

0049  (1.24) 

0  050 

(1  25) 

Slow  1)  perf  (v.) 

0  767  (103) 

0  620 

(104) 

Fat  (>>) 

0  049  (125) 

0  230 

(109) 

Cardiac  output  (1  It)  (^,) 

M3‘  (108) 

348' 

(1  12) 

Waling  value 

- 

486' 

(1  12) 

Alveolar  ventilation  -ate  (1  Ti) 

164  (111) 

288' 

(I  50) 

Waling  value<^^) 

- 

683' 

(126) 

Comparlmenl  perfusions  0.'’b) 

Liver  (<j,) 

0.282  ('  2<) 

906 

(135) 

RapidI)  perf  (^.1 

0.S76  (124) 

153 

(125) 

Slow!)  perf  (<? ' 

0.120  (101) 

87.0 

0  04) 

Waling  value 

- 

225 

(I  17) 

F*'  (lA 

0  J02  (125) 

17.4 

0  09) 

Partition  coefficients  (unilless) 

Blood  gas  LptI 

169  (1.97) 

110 

0  97) 

Liver.'blood  (p,) 

301  (269) 

505 

(9.37) 

Liver  gas  (p^) 

50.9  (1.97) 

60.6 

(*3«) 

Liver  gas  (p„) 

50.9  (1.97) 

60.6 

(8  36) 

Rapid  "blood  (p.) 

3.01  (4.14) 

5  05 

(5.69) 

Rapid  gas  (/I,) 

50.9  (3  51) 

606 

(4.92) 

Slow.'blcxxl  ip,) 

2  59  (2.54) 

2  66 

(110) 

Slow  gas  (p,) 

43  J  (1.97) 

319 

(101) 

Fat, "blood  ipf) 

4*  3  (2  56) 

102 

(2.15) 

Fil/gas  ipj,) 

816  (1.93) 

1230 

(2.15) 

Metabolic  constants 

F— -  (mg/b) 

3%6  (283) 

0J3 

(2.84) 

K,  (mg/1) 

1.47  (124) 

1.86  (12.3) 

'The  uncertainly  factors  are  estinuted  frocn  1000  LaliD-bypcTcvbe  samples  such  thal  for  preferred  value 
the  inlervaJ  from  A^UF  to  A,  x  UF  coouitu  95%  of  the  values. 

‘Preferred  values  for  mtoe.  2a  h/d. 

'Preferred  values  for  sleeping  buniaus,  S  b/d. 

'Preferred  values  for  waking  humans,  16  b/d. 
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al  ihe  definition  of  mode!  parameters  and  the  units  adopted  follow  Reitz  and  Nolan 
(51  Absolute  compartment  solumes  (»,)(in  liters)  are  related  to  compartment  propor¬ 
tions  Table  I)  by, 

r,=  *  bs\  (i  =  l,r.sj) 

uere  is  body  \scight  in  kilograms.  The  maximal  metabolic  rale  (k'mn)  is  scaled 
according  to  k'm,,  =  k',r.jw  (5)  The  and  parameters  are  those  for 

which  uncertainties  are  estimated. 

The  partition  coefficients  that  are  measured  directly  in  vitro  are  limited  to  some 
of  the  tissue, 'gas  coefficients  {pt  and  p,,.  i=l,rjJ).  Tissue/hlood  partition  coefficients, 
which  are  those  actually  used  in  the  PBPK  model,  are  estimated  by  dividing  the 
tissue  gas  coefficients  by  the  blood  'gas  coefficient  {p,  =  PiJpb,  i  = 

The  dose  surrogates  that  are  estimated  herein,  and  for  which  uncertainties  are  esti¬ 
mated,  arc  (I)  the  avciagt  daily  area  uiidcr  the  liver  concentration-time  curve  for 
the  parent  (AUCL).  (2)  the  average  daily  area  under  the  arterial  blood  concentra¬ 
tion-time  curve  for  the  parent  (AUCA).  and  (3)  the  average  daily  amount  of  parent 
metabolized  per  volume  of  liver  tissue  (CML). 

PARAMETER  UNCERTAINTY  EVALUATION  FOR  THE  TETRACHLOROETHYLENE 
MODEL 

For  each  parameter  we  have  identified  a  preferred  value  (Table  I),  and  have  speci¬ 
fied  a  probability  distribution  to  represent  uncertainty.  The  preferred  value  repre¬ 
sents  some  summary  of  information  available  from  the  literature  (see  below).  The 
probability  distribution  is  selected  in  such  a  way  as  to  assign  relatively  high  probabili¬ 
ty  to  values  that  are  close  to  the  preferred  value.  A  useful  device  for  communicating 
uncertainties  is  to  define  an  ‘uncertainty  factor’  (UF  ?  1)  such  that  for  a  preferred 
parameter  value  {hp)  the  range  of  values  /ip  TJF  to  hp  x  UF  is  considered  highly  prob¬ 
able  (we  assume  probability  95%).  Table  I  gives  uncertainty  factors  estimated  from 
the  observed  distributions  of  the  values  generated  for  the  parameters,  which  are 
based  on  the  probability  distributions  desenbed  below. 

Specification  of  probability  distributions  to  express  uncertainty  in  parameter  val¬ 
ues  involves,  first  of  all.  selection  of  appropriate  families  of  distributions  for  individu¬ 
al  parameters,  or  in  some  cases  for  groups  of  inierrelalcd  parameters.  The  families 
selected  for  expressing  uncertainty  are  the  log-normal  or  truncated  log-normal  family 
(for  partition  coefficients  and  metabolic  constants)  and  the  Dirichlet  family  (for  com¬ 
partment  volumes  as  proportions  of  total  body  volume).  Other  PBPK  model  param¬ 
eters  are  functionally  related  to  the  compartment  volumes  and  their  empirical  distri¬ 
butions  depend  on  the  distribution  of  the  volumes. 

Uncertainty  factors  wcie  derived  largely  on  the  basis  of  analyses  of  variation  in 
reported  measurements  or  estimates.  That  being  the  case,  different  sources  of  uncer¬ 
tainty  and  variation  must  be  recognized.  Variation  among  average  measurements  re- 
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poncd  in  difTcfcnt  cxpcnmcnts  represents  a  different  source  or  sources  of  uncerlainl) 
than  sarialion  from  measurements  tai-cn  on  different  individuals  in  the  same  experi¬ 
ment.  Furthermore,  \shen  a  parameter  salue  is  not  measured  directly  so  that  a  value 
of  another  variable  is  attributed  to  that  parameter,  additional  uncertainty  is  in¬ 
troduced  Consequently,  the  evaluation  of  total  uncertainty  theoretically  requires 
combining  distributions  represer.tirtg  multiple  levels  of  uncertainty. 

Of  interest  for  test  species  PBPK  parameters  is  uncertainty  with  respect  to  an  aver¬ 
age  parameter  value.  This  is  the  case  because  an  aggregate  response  variable  (propor¬ 
tion  of  animals  with  tumor)  is  to  be  related  to  an  aggregate  predictor  (the  single  dose 
surrogate  value  assumed  for  a  treatment  group).  For  humans,  in  contrast,  interest 
is  in  individual  variation  of  dose  surrogate  levels  dve  to  variation  in  PBPK  model 
parameter  values.  Therefore,  for  humans,  it  is  desirable  to  incorporate  into  an  esti¬ 
mate  of  relevant  total  uncertainty  a  component  representing  inter-individual  varia¬ 
tion  specifically.  Consequently  the  probability  distributions  derived  below  for  hu¬ 
mans  appl,.  to  individuals  selected  at  random,  but  distributions  derived  for  rodents 
apply  Vo  mean  s  alucs  for  groups  of  animals  in  the  bioassay  considered. 

Presented  here  are  unceriointy  derivations  for  mice  and  humans.  The  concepts  ex¬ 
emplified  by  the  analysis  for  mice  generalize  to  other  test  Sijecies. 

Partition  coefficients 

Preferred  values  for  tissuc/gar  partition  coefficients  (equated  with  the  medians  of 
the  log-normal  distributions  charactcrio'ng  uncertainties)  are  taken  from  Reiut  and 
Nolan  [5]  The  blood  gas  preferred  value  (12)  is  the  geometric  mean  of  the  value  from 
Reitz  and  N  Ian  [5]  (10  3)  and  3  valuer  (9  1,  13.1,  18.9)  rejxirtcd  by  Fiserova-Berge- 
rova  [6]. 

Three  sources  of  uncertainty  arc  recognized,  one  or  more  of  which  arc  relevant 
to  the  esii. nation  of  total  uncertainty  for  a  given  partition  coefficient.  For  parameters 
measured  directly  by  vial  equilibration  techniques,  uncertainty  is  due  to  differences 
in  values  estimated  in  different  experiments  (source  1).  These  parameters  are  the  ro¬ 
dent  blood  gas,  liver/gas,  fat, 'gas,  and  slow  gas  (actually  musclc/gas)  and  the  human 
blood,  gas  coefficients.  The  unmeasured  paiamcters  -  the  rapid 'gas  coefficient  in  ro¬ 
dents  (equaled  to  the  rodent  liver  gas  coeffiaent)  and  all  tissue,  gas  coefficients  other 
than  blood  gas  in  humans  (equated  to  the  average  of  the  corresponding  coefficients 
estimated  in  rats  and  mice)  -  have  additional  uncertainty  due  to  the  attribution  of 
the  measured  values  to  those  parameter  values  (source  2).  For  all  human  coefficients, 
individual  variation  is  also  relevant  (source  3). 

Each  source  of  uncertainly  is  represented  by  a  separate  log-normal  distribution 
with  specified  geometric  standard  deviation  (GSD).  Therefore,  the  distribution  de¬ 
scribing  total  uncertainty  is  also  log-normal,  with  geometric  variance  (squared  geo¬ 
metric  standard  deviation)  given  by  the  sum  of  geometric  variances  expressing  specif¬ 
ic  contributing  uncertainties.  Details  of  the  derivation  of  source-specific  uncertainty 
distributions  arc  available  from  the  authors. 


VI-1-92 


ll  is  commor  practice  in  unceriainiy  analyses  tc  dtrise  <?i5lribulion  parameters  by 
directly  assigning  subjective  probabilities  to  ranges  of  parameter  values  |7].  An  analo¬ 
gous  procedure  has  been  used  here  to  modify  the  uncertainty  estimates  for  the  parti¬ 
tion  coefficients. 

Plausible  constraints  on  values  of  partition  coefficients  are 
p,oU~  'shere  />,.  and  p,o  are  the  water/gas  and  oil  gas  partition  coefficients, 

respectively.  Using  the  preferred  values  (Table  I)  and  source-specific  uncertainty  fac¬ 
tors  as  discussed  above  to  define  the  log-normal  distributions,  the  value  of  p^j  will 
be  larger  than  the  value  of  1917  estimated  for  p^o  by  Sato  and  Nakajima  (8),  with 
probably  about  20%.  Similarly,  in  sampling  from  log-normal  distributions  for  the 
partition  coefficients,  values  of  other  tissue/gas  partition  coefficients  were  occasional¬ 
ly  larger  than  values  of  the  fat/gas  coefficient.  (Other  suggested  constraints  were  vio¬ 
lated  with  very  low  probability.)  Consequently,  the  distribution  of  the  fat/gas  parti¬ 
tion  coefficient  was  truncated  at  1917,  and  other  tissue,  gas  coefficients  were  truncated 
at  the  value  of  the  fat/gas  coefficient.  These  adjustments  represent  a  significant  altera¬ 
tion  of  both  the  shape  and  the  variance  of  the  distribution  representing  uncertainty 
of  the  fat,  gas  partition  coefficient,  and  relatively  slight  alteration  of  the  distributions 
for  other  tissue,'g3s  partitior'  coefficients. 

Compartment  volumes 

For  each  species,  uncertainty  regarding  compart.ment  proportions  of  total  body 
volume  is  expressed  using  the  Dirichlet  distribution.  For  a  set  of  random  complemen¬ 
tary  proportions,  the  Dirichlet  distribution  function  can  be  expressed  in  terms  of  the 
expected  proportions  and  a  parameter  6  which  determines  the  variances  and  co-var- 
iaiiCvS  of  the  ra  i  om  proportions  (9).  Consequently,  it  is  convenient  to  equate  the 
preferred  proportions  v^^  (i  =  l,rjJ)  given  in  Table  I  to  the  expectations  of  the  cone- 
sponding  random  proportions  =  The  value  chosen  for  9  expresses  the 

uncertainty  regarding  the  joint  distribution  of  the  proportions;  more  specifically  var¬ 
iance  of  v^  is  given  by; 

Var(v  J  =  (1  -  (0  -I-  I ),  f  =  (I ) 

Fixing  the  variance  of  one  compartment  proportion  determines  9.  Variation  in 
liver  volume  had  been  used  because  prebminary  analysis  revealed  that,  for  the  3  dose 
surrogates  studied  and  for  inhalation  exposure,  the  model  is  roughly  equally  sensitive 
to  all  compartment  volumes  (for  the  AUCL  and  AUCA  surrogates)  or  more  sensitive 
to  the  liver  volume  than  to  other  volumes  (for  the  CML  surrogate).  Also,  pubUshed 
information  is  easily  incorporated  since  the  compartment  is  identified  unambiguously 
with  a  specific  organ. 

On  the  basis  of  bver  volume  measurements  for  mk*  reported  in  Arms  and  Travis 
(101,  the  preferred  value  of  tif,  and  9  for  mice  arc  estimated  to  be  0.0S6  and  1456, 
respectively.  Volume  proportions  for  other  compartments  are  taken  from  Reitz  and 
Nolan  (S],  adjusted  for  the  change  in  v 
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For  humans,  ihc  paramcicr  value  of  »v  (0  0?6)  is  also  taken  from  Arms  and  Tra' 
[10]  and  other  volume  proportions  from  ReiU  and  Nolan  |5],  adjusted,  again,  for  t 
diflTerencr  bct\vccn  vj,  values  from  the  tv^o  sources.  For  0,  individual  variation  shou 
be  considered  (as  discussed  above).  Data  on  individual  variation  in  liver  volume 
humans  (as  a  proportion  of  body  weight)  is  not  available.  Caster  et  al.  [11]  repc 
variances  of  vj,  tiiimates  for  individual  rats,  which  we  assume  to  express  plausit 
levels  of  individual  variation  for  humans.  Using  value  from  Caster  et  al.  [1 1]  in  Equ 
tion  (1),  the  solution  for  S  is  1621.  We  tentatively  adopt  this  value  to  represent  unce 
tainty  for  humans,  noting,  however,  that  individual  differences  do  not  account  f 
total  uncertainty  regarding  compartment  volumes,  so  that  some  under-estimate 
uncertainty  is  possible. 

Circulatory  and  \ertlilalory parameters 

The  parameters  <^e<^r'''g  the  ».oriipartrncnt  proportions  of  body  volume,  cardi 
output  (i^.),  ventilatory  output  and  perfusion  rales  of  specific  compartments  (i 
i-l.r.s/)  are  functionally  interrelated.  An  account  presented  here  demonstrates  t! 
sort  of  physiological  considerations  that  are  relevant,  along  with  possibilities  f 
quantifying  the  way  that  these  considerations  may  determine  uncertainty.  Data  r 
garding  circulatory  and  ventilatory  parameters  for  rodents  and  humans  have  bci 
summarised  by  Arms  and  Travis  [10]. 

Distributions  for  qp,  and  the  q,'s  are  derived  as  follows: 

(1)  Given  a  value  sampled  from  the  Dirichlet  distribution  and  the  correspondii 
absolute  compartment  volume  (vj,  q,  (the  absolute  compartment  perfusion  rate) 
given  by  q,  =  (v,  x  q^p  x  qtf)/v,p  where  is  the  preferred  proportion  of  total  cardi; 
output  ''  cctcd  to  compartment  I'llO]  and  qtp  is  the  preferred  total  cardiac  outpi 
TThc  total  cardiac  output  is  then  given  by  =  ?/+?, +  ?,  +  ?/.  These  values  of  lot 
and  compann.ent-spev;ific  perfusion  are  assumed  to  hold  for  sleeping  humans  at 
for  rodents  during  all  hours. 

(2)  For  humans,  sleeping  and  waking  values  of  cardiac  output  and  alveolar  venlil 
lion  rate  are  desired  The  ‘waking’  value  {qc»)  is  computed  by  qn  =  qr  +  {qne- 
q^f)  X  edq^  where  q,  is  the  'sleeping'  value  computed  in  step  1,  q^p  and  q^r  are  pre¬ 
ferred  values  for  sleeping  and  w  aking  individuals,  and  edq^  is  distributed  log-normal¬ 
ly  with  preferred  value  1  and  uncertainly  factor  lobe  specified.  The  increment  to  total 
cardiac  volume  is  assumed  to  be  directed  entirely  to  the  slowly  perfused  compart¬ 
ment.  PrcferTcd  values  of  waking  and  sleeping  cardiac  output  and  alveolar  ventila¬ 
tion  rate  are  derived  from  Hattis  et  al.  W- 

An  uncertainty  factor  for  edq,  is  also  based  on  the  data  presented  in  Hattis  et  al. 
[4].  An  arbitrarily  low  probability  (0.001)  is  assigned  to  an  increment  in  q,  values 
(qcm~R<)  as  large  as  the  difference  between  values  for  sleeping  and  light  exercise  (234 
1/h).  Using  this  difference,  the  assumed  difference  in  preferred  values  {q„p- 
q^P  -  1 38),  the  log  normal  distribution  assumed  for  edq^,  and  the  relationship  be- 
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IwccT  ihe  GSD  of  ihc  dislribulion  and  ihc  associalfd  UF,  ihc  uncertainly  factor  of 
13  is  obtained. 

(3)  Ventilatory  volume  (q^)  is  assumed  to  be  related  linearly  to  the  cardiac  volume 
(^,).  For  humans.  -  <7;,^  =  [(<?;,.  a -<?;./-) -?./•)])»  (^r  -  q^r)  The  same  relation¬ 
ship  holds  for  waVing  values,  <7^.  and  q^^  For  mice,  an  analogous  treatment  requires 
specification  of  ?  second  <7.-?^  coordinate  (to  replace  q,.f  and  qp,f),  in  addition  to 
the  values  representing  ordinary  activity  levels  Data  from  anesthetized  mice  have 
been  used  to  determine  the  second  coordinate  (lOj. 

Metabolic  constants 

After  scaling  estimates  (to  get  values)  and  express  ng  in  comparable 
units,  geometric  means  of  the  values  reported  in  Hattis  et  al.  (4)  and  Reitz  and  Nolan 
(5]  arc  equaled  with  the  preferred  values 

It  is  assumed  here  that  the  uncertainty  regarding  the  metabolic  parameters  arises 
chiefly  from  different  exjxrimental  approaches  and  assumptions  applied  in  deriving 
estimates.  It  may  be  argued  that  an  uncertainty  factor  based  on  a  geometric  standard 
error  is  an  appropriate  measure  of  uncertainty,  the  geometric  standard  error  is  relat¬ 
ed  inversely  to  the  sample  size,  representing  formally  the  inverse  relationship  between 
uncertainty  and  additional  information.  However,  such  a  formalization  is  based  on 
the  assumption  that  the  available  measurements  are  independent  and  identically  dis¬ 
tributed  (i.i  d.).  The  estimates  of  metabolic  constants  are  suspected  to  deviate  strong¬ 
ly  from  these  assumptions.  A  likely  pattern  of  deviation  from  the  i.i.d.  case  is  for 
clusters  of  estimates  based  on  similar  procedures  to  assume  similar  values.  Such  a 
pattern  can  lead  to  too-small  estimates  of  uncertainty.  Here,  some  conservatism  (im¬ 
plying  greater  uncei  .dinty)  is  introduced  by  basing  an  uncertainty  factor  on  the  geo- 
metne  s.a.itja.'d  .ether  tttar*  tiic  geometric  stdiiuaro  ciiOi.  In  other  wv>ios, 

uncertainly  regarding  the  metabolic  constants  is  represented  by  a  Icfc-normal  distri¬ 
bution  having  the  same  geometric  standard  deviation  as  the  sample  geometric  stand¬ 
ard  deviation  computed  from  the  av  ailable  estimates  (4,  5]. 

A  tentative  uncertainty  factor  for  is  given  by  exp  (2  x  0.532)  =  2.9,  for  mice 
and  humans,  where  0.532  is  z  GSD  denved  as  a  weighted  average  of  species-specific 
GSD  values.  Uncertainty  for  is  evaluated  in  a  similar  manner,  resulting  in  an  un¬ 
certainty  factor  of  13  for  mice  and  humans  The  representation  of  uncertainty  with 
respect  to  these  metabolic  constants  is  considered  problematic,  especially  because 
their  estimation  via  optimization  of  the  model  fit  to  in  vivo  data  is  dependent  on  esti¬ 
mates  for  the  other  model  parameters  The  ‘Discussion’  (sec  below)  elaborates  on 
this  point. 

PROPAGATION  OF  PARAMETER  UNCERTAINTIES 
Procedure 

Parameter  distributions  defined  in  the  section  on  ’Parameter  UrKcriainty  Evalua- 
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lion'  for  mice  and  humans  uere  sampled  using  ihe  Lalin-h>"percube  melhod  |7,  12). 
For  each  random  selection  of  parameter  salues,  the  operations  depicted  in  Figure 
2  skcre  performed  to  generate  distnbutions  of  dose  surrogate  values  and  risk  esu- 
mates. 

In  order  to  charactenze  uncertainty  in  the  3  dose  surrogates  (AUCL,  AUCA, 
CML),  100  sets  of  parameters  were  generated  and  the  PBPK  model  was  run  for  each 
set.  Mice  were  exposed  as  in  the  carcinogenicity  bicassay  via  inhalation  at  0,  100 
and  200  ppm,  for  6  hours  per  day  in  5  consecutive  days  per  week  (2).  Humans  were 
assumed  to  be  exposed  to  50  ppm  (the  current  OSHA  standard)  for  8  hours  per  day 
in  5  consecutive  days  per  week. 

Uncertainty  in  dose  surrogate  values  was  propagated  further  to  evaluate  uncer¬ 
tainty  in  human  risk  estimates.  A  single  tumor  response,  hepatocellular  carcinoma, 
is  considered  here  for  illustrative  purposes. 
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Surrogate  dose  values  were  related  to  the  tumor  response  rales  using  a  version  of 
oLOB  M82  [13]  to  implement  the  multiviagc  model.  For  the  estimation  of  risk  un¬ 
certainty,  tumor  response  rates  were  also  considered  uncertain.  For  each  of  100  ran¬ 
domly  generated  sets  of  dose  surrogate  values  (one  dose  surrogate  value  per  treat¬ 
ment  group)  the  multistage  nvodel  with  a  fixed  number  of  stages  was  refitted,  based 
on  those  dose  surrogate  values  and  random  tumor  counts  for  the  treatment  groups 
generated  from  a  binomial  distribution  with  parameters  S,  and  p,.  Here  S,  is  the 
nurriLicr  of  animals  in  treatment  group  r  (48  at  0  ppm,  50  at  100  and  2(X)  ppm),  while 
p,  IS  the  estimated  response  proportion  for  dose  group  /  obtained  by  fitting  the  multi¬ 
stage  model  to  the  observed  response  rates  (I  at  0  ppm,  13  at  lOO  ppm,  and  36  at 
200  ppm)  and  the  dose  surrogate  values  corresp'onding  to  the  preferred  mouse  pa¬ 
rameter  values. 

Finally,  a  distribution  of  human  maximum  liV-clihood  and  upper  bound  estimates 
of  riskt  were  generated  by  pairing  each  of  the  KX)  sets  of  multistage  model  parameters 
with  a  randomly  selected  human  dose  surrogate  value  corresponding  to  the  human 
exposure  scenario  of  interest  Extra  risk  is  defined  by  /?  =  (/*(d)  -  A0))/(1  “  AO)), 
when  P{<f)  IS  the  lifetime  probability  of  observing  a  tumor  given  exposure  that  results 
in  a  surrogate  dose  value  of  d.  No  adjustments  were  made  to  the  risk  estimates  ob¬ 
tained  in  this  nicir  ncr  Thus,  for  example,  no  adjustment  was  made  to  account  for 
the  somewhat  different  proportions  of  the  human  and  mouse  lifespans  lived  prior  to 
first  exposure. 

RESULTS 

Table  II  gives  se.^cted  percentiles  for  the  distnbulions  of  human  dose  surrogate 
values  and  extra  nsks  Median  surrogate  values  are  approximately  equal  to  the  dose 
surrogate  values  computed  using  preferred  parameter  values,  which  are  23.9  for 
AUCA.  1 18  for  AUCL.and  27  8  for  CML. 

Median  nsk  estimates  vary  substantially  among  the  3  surrogate  doses.  For  com¬ 
parison,  consider  the  nsks  estimated  in  the  ’traditional’  manner,  using  applied  doses 
with  no  pharmacokinetic  transformations  (Table  III)  Maximum-likelihood  risk  esti¬ 
mates  (MLE)  obtained  without  using  a  PBPK  model  (assuming  mice  and  humans 
are  equally  susceptible  to  doses  expressed  in  mg.'kg  body  wt.  per  day)  are  higher  than 
maximum  likelihood  estimates  based  on  metabolite  in  bver,  but  lower  than  MLEs 
based  on  parent  concentrations  For  tctrachloroethylene,  it  appears  that  the  structur¬ 
al  uncertainty  associated  with  the  selection  of  an  appropriate  dose  metric  for  cross- 
species  extrapolation  is  of  relatively  greater  importance  than  is  the  uncertainty  asso¬ 
ciated  with  the  values  of  the  PBPK  model  parameters. 

Also  shown  in  Table  III  are  the  risk  estimates  ob'iained  by  using  the  preferred 
mouse  and  human  PBPK  parameters  to  estimate  risk  (i  *•.  the  risk  estimates  that 
would  be  obtained  if  the  parameter  values  were  known  without  uncertainty  to  be 
equal  to  the  preferred  values).  Except  for  the  estimates  based  on  CML,  the  median 
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TABLE  il 

SLLEC  TLD  PL  R<.  ENTILES  FROM  THE  SIMULATED  DISTRIBUTIONS  OF  HUMAN  DOSE 
SURROGATE  AND  RISK  VALUES  (BASED  OS  100  LATIN  HYPERCUBE  SAMPLES) 

Soticfilt  Prtctnlilei 

25»  50*  97  5* 


fllMC  Su"OgillfS 

AUCA 

dai!)  jica  undfr  (he  arlenal 
c-jr'f,  mg  h' '  1' ') 

aucl 

(a'-ftage  dail>  J'ca  undti  iht  li^ef  con- 

crn'.talic-n  cjn-c.  mg  h'  'T' ') 

UML 

(a'-<--.i('^  d.iil)  amc  '^n;  me (aS.'-li'cd  per 
liCt  ^oGme,  mg  II 

Hji.p -m  !:ii\',hr<od  r!:i.'rijies  '  MLLsi  Jot 
fii'a  'ui" 

AUCA 

aucl 

CML 

•E  ora  mk  /TO))  (I  -  /TO)).  »herc  /\(f)  is  the  bfeUme  probabdily  of  tumor  when  npo«d  to  do* 

J  Here  dee  rrnponds  lo  ar  8  h  d.  5  d  »  expoturt  to  50  ppm 

table  ill 

RISK  f  STiMATES*  obtained  Vk  ITHOUT  USE  OF  A  PBPK  MODEL,  AND  USING  THE  PBPK 
MODEL  'AiTH  NO  U.  '-ERT AINT't' 


;  t 

18  9 

24  4 

303 

43  8 

3  1 

54  9 

116 

243 

1060 

2  24 

13  9 

27  8 

St.8 

117 

0  0275 

0  0910 

0  164 

0  243 

0  494 

^391 

0  244 

0  480 

0  700 

0  982 

r65E-9 

3  73E-5 

8  50E-4 

7  06E-3 

0  049' 

Anal>us 


Risk  nlimaics 


MLE  Llppcr  bound 


No  PBPR  '.rantforma'.iOO* 

PBPK  iranxfonT.aiior  '*ilh  no  unccruinly 

AUCA 

AUCL 

CML 


5  57  *  10-' 

0  126 
0d23 

1.95*  10- » 


4  28  *  10-> 

0  238 
0506 

7.00*  10-' 


•E»lra  niks  for  8  h/d.  5  d.**  exposure!  lo  50  pptn. 

‘Humans  and  mset  are  assumed  to  be  equally  susceptible  lo  telrachloroethyltoe  when  administered  does 
art  expressed  in  mg  kj"  '-d* 


of  ihc  simulated  risk  estimates  is  dose  lo  ihe  risk  estimale  obtained  with  no  uncer- 
lainty  In  the  case  of  CML,  risks  terxJ  lo  be  skewed  to  the  right  in  comparison  with 
the  estimate  corresponding  to  the  preferred  parameter  values.  This  is  a  result  of  the 
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fact  that  the  linear  term  of  the  multistage  model  fit  to  the  preferred-value  dose  and 
response  data  is  zero.  Alteration  of  CML  doses  and  response  rales  introduced  in  the 
course  of  the  uncertainty  analysis  cannot  decrease  the  linear  term  and  this  would  tend 
to  increase  the  iisk  estimates  at  the  relatively  low  dose  corresponding  to  the  human 
exposure  scenano.  Nk'ilh  the  AUCL  and  AUCA  surrogates,  the  linear  terms  are  posi¬ 
tive  when  the  multistage  model  is  fit  to  the  observed  response  rates  and  doses  corre¬ 
sponding  to  the  preferred  parameter  values. 

DISCUSSION 

Evaluations  of  the  carcinogenicity  of  tetrachlorocthylene  have  implicated  a  meta¬ 
bolite  as  the  moiety  responsible  for  the  tumorigenic  responses  following  exposure 
(14).  The  high  degree  of  non-linearity  between  the  CML  dose  surrogate  and  the  hepa¬ 
tocellular  carcinoma  response  rates  in  female  mice  (as  discussed  in  the  preceding  sec¬ 
tion)  is  interesting  in  that  light.  Either  there  are  non-linear  steps,  in  addition  to  the 
saturable  metabolic  transformation,  leading  to  the  interaction  that  initiates  the  carci¬ 
nomas  (i  e  .  the  surrogate  represented  by  CML  is  not  close  enough  to  the  events  caus¬ 
ing  cancer  to  yield  linearity  between  that  measure  of  dose  and  the  response),  or  the 
PBPK  model  is  not  adequate  for  describing  tetrachlffroethylene  pharmacokinetics. 
We  are  continuing  investigation  of  tetrachlorocthylene  cardnogenicity  in  hopes  of 
eluddating  the  substance  and  the  mechanisms  responsible.  This  may  entail  the  devel¬ 
opment  of  alternative  means  of  characterizing  and  calculating  the  risks. 

The  approach  to  uncertainty  analysis  that  has  been  illustrated  in  the  context  of 
tetrachlorocthylene  risk  assessment  has  several  features  that  are  essentia!  for  an  ap¬ 
propriate  uncertainty  assessment.  First,  the  uncertainties  related  to  the  PBPK  model 
parameters  are  formally  expressed  in  terms  of  probability  distributions  that  can  in¬ 
corporate  multiple  levels  (or  sources)  of  uncertainty.  The  parameter  uncertainty  dis¬ 
tributions  are  based  on  the  inspection  and  analysis  of  relevant  data,  sometimes  in 
combination  w  ith  the  application  of  more  subjective  evaluations  (or  those  based  on 
expert  opinion)  of  reasonable  bounds  for  particular  parameter  values.  Second,  the 
parameters  are  not,  in  general,  regarded  as  being  independent  of  one  another.  The 
implied  correlation  structures  (as,  for  example,  among  the  volumes  of  the  compart¬ 
ments)  arc  modeled  using  multivanate  distributions.  Moreover,  certain  assumed 
functi.  nal  relationships  between  the  parameters  have  been  maintained,  even  though 
the  magnitudes  of  those  relationships  (the  ratio  of  cardiac  and  pulmonary  rates,  for 
example)  are  subject  to  ui’ceriainty.  Finally,  the  approach  employs  an  efficient  tech¬ 
nique  (Monte  Carlo  stimulation  using  the  Latin-hypcrcube  procedure)  to  make  ex¬ 
plicit  the  relationships  between  parameter  uncertainties  and  uncertainties  in  dose  sur¬ 
rogate  values  and  risk  estimates. 

By  far  the  most  difficult  aspect  of  this  approach  is  the  characterization  of  the  un¬ 
certainties  in  the  parameters.  This  requires  extensive  review  of  the  literature  and  at¬ 
tention  to  the  different  sources  of  uncertainty  that  are  inherent  in  the  observed  varia- 
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lions  in  measurements.  It  is  appropriate  to  stress  that  the  estimates  of  uncertainty 
derived  here  are  tentative.  First  of  all.  as  illustrated  by  the  measurement  of  individual 
variation  in  liver  volume  as  a  proportion  of  body  weight,  various  of  the  GSD  csb- 
mates  arc  not  based  on  data  that  are  directly  relevant,  but  on  data  from  contests 
that  are  considered  more  or  less  analogous  to  the  context  of  immediate  interest.  Mo¬ 
reover,  the  variations  identified  with  specific  sources  of  uncertainty  may  not  reflect 
only  those  sources  of  uncertainty.  For  example,  individual  variation  will  contribute 
to  observed  variation  between  mean  measurements  reported  by  different  laborato¬ 
ries,  particularly  if  the  number  of  animals  per  mean  is  small. 

Certain  parameters  are  difficult  to  measure  and  equally  difficult  to  charactenze 
with  respect  to  their  uncertainty  .  WTicn  parameter  estimates  are  based  upon  data  col¬ 
lected  in  vivo,  as  is  the  case  with  metabolic  constants,  the  values  obtained  are  deter¬ 
mined  in  pan  by  the  values  assumed  for  other  physiological  parameters,  as  well  as 
by  the  structural  assumptions  underlying  the  PBPK  model.  Also,  uncertainties  in  es¬ 
timates  of  parameters  fitted  simultaneously  by  an  optimization  program  are  not  inde¬ 
pendent,  because  the  fit  of  the  model  to  the  data  is  determined  jointly  by  all  of  the 
fitted  parameters.  Information  on  the  correlation  among  estimates  of  such  param¬ 
eters  is  necessary  for  a  complete  treatment  of  their  uncertainty.  This  information  is 
not  currently  available,  and  so,  for  example,  and  arc  treated  as  though  they 
were  independently  determined.  Because  of  these  difficulties,  the  evaluation  of  uncer¬ 
tainly  for  metabolic  constants  could  be  improved  by  techniques  that  account  for  the 
uncertainty  in  the  other  parameters  and  that  simultaneously  estimated  uncertainty 
for  ^mxx  and  K^n- 

Other  icfinemcnts  to  the  characterization  of  uncertainties  in  PBPK  model  param¬ 
eters  are  p^^siblc.  In  particuhr,  the  use  of  variance  components  analysis  to  denve 
more  pure  estimates  of  source-specific  uncertainties  is  of  interest. 

The  results  of  the  analysis,  the  distribution  of  surrogate  doses  and  of  risk  estimates, 
have  several  potential  uses.  The  distribution  of  simulated  risk  estimates  can  be  inter- 
preled  as  follows,  keeping  in  mind  the  fact  that,  precisely  because  of  uncertainty  with 
respect  to  the  parameter  values  and  individual  variation,  consideration  of  a  single 
nsk  estimate  is  not  adequate.  The  proportion  (F)  of  the  simulated  estimates  that  fall 
in  an  interval  (7)  may  be  interpreted  as  the  probability  that  the  true  risk  for  a  random¬ 
ly  selected  individual  is  in  /,  when  the  uncertainties  are  taken  into  account.  This  inter¬ 
pretation  should  be  useful  in  risk  management  decision-making  contexts.  One  major 
advantage  of  this  type  of  analysis  0  risk  assessment  that  uses  PBPK  modeling  and 
considers  input  uncertainties)  is  that  reasonable  variations  in  risk  estimates  become 
explicit.  Traditionally,  using  administered  doses,  the  uncertainties  associated  with  the 
point  estimates  of  risk  have  not  been  emphasized. 

The  distribution  of  the  dose  surrogates  and  the  risk  estimates  result  from  uncer¬ 
tainties  concerning  true  (average)  values  as  well  as  variation  around  those  averages. 
The  uncertainty  regarding  the  true  values  is  theoretically  reducible  via  further  experi¬ 
mentation.  The  approach  illustrated  here  can  assist  in  directing  the  efforts  to  reduce 
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uncertainty  If.  for  example,  a  set  of  parameters  is  fixed  at  an  assumed  axerape  value, 
the  reduction  in  the  spread  of  the  risk  estimates,  hclow  (hat  obtained  when  all  param¬ 
eters  arc  allowed  to  vary,  provides  an  indication  of  the  potential  value  of  additional 
information  with  respect  to  that  set  of  parameters  Such  analyses  lake  into  account 
both  sensitivity  of  model  output  to  the  values  assumied  for  the  parameters,  and  the 
current  level  of  uncertainty  regarding  the  parameters 

The  approach  to  uncertainty  analysis  presented  here  has  been  illustrated  for  a  sim¬ 
ple  PBPK  model  for  letrachloroeihylenc  It  is  applicable  to  PBPK  models  of  greater 
complexity  as  well.  Some  of  the  parameters  of  those  models  -  such  as  partition  coeffi¬ 
cients,  compartment  volumes,  and  flow  rates  -  are  identical  to  those  in  the  tetrachlo- 
rocihylenc  model  used  herein  Thus,  for  those  parameters,  the  ground  work  has  been 
laid  here  for  their  analysis  In  genera),  the  considerations  described  here  are  mean¬ 
ingful  and  useful  for  any  analysis  of  PBPK  model-based  nsk  assessment. 
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VOLUME  VI 


SENSn  IVITY/UNCERTAINTY  ANALYSIS 

PART  2  OF  2  PARTS 

CONTRIBUTIONS  OF  INDIVIDUAL 
PARAMETERS  TO  OVERALL  UNCERTAINTY 


A.  INTRODUCTION 


In  Part  1  of  this  volume,  a  sensitivity  analysis  for  individual 
parameters  of  a  simple  PERC  PBPK  model  was  completed.  Appendix  VI-l-A 
presented  an  uncertainty  analysis  for  chat  same  model.  The  uncertainty 
analysis  built  upon  the  discussion  in  Part  1  concerning  parameter 
uncertainties.  The  uncertainty  analysis  also  presented  results  in  terms  of 
distributions  of  output  variables,  in  particular  distributions  of  risk 
est imates . 

In  this  part  of  the  volume,  one  other  sec  of  results  is  presented. 

Those  results  also  apply  to  the  uncertainty  analysis  and  again  focus  primarily 

on  the  distribution  of  output  variables.  In  this  case,  however,  the  analysis 
highlights  the  contributions  of  individual  parameters  or  sets  of  interrelated 
parame  ters . 

The  goal  of  the  analysis  reported  here  was  to  determine  which  parameters 
contributed  most  significantly  to  the  variability  in  the  dose  surrogate  or 
risk  estimates.  A  variable  can  contribute  to  output  uncertainty  through  a 
combination  of  1)  the  sensitivity  of  the  models  to  the  output,  and  2)  the 
degree  of  uncertainty  (variability)  associated  with  the  parameter  estimate.  A 
model  is  said  to  be  sensitive  to  a  parameter  if  relatively  small  changes  in 

the  parameter  values  result  in  relatively  large  changes  in  model  output 

values.  As  stated  in  Section  1  of  Part  1  of  this  volume,  a  parameter  with  a 
high  degree  of  uncertainty  associated  with  it  does  not  contribute  greatly  to 
output  uncertainty  if  the  model  is  relatively  insensitive  to  that  parameter. 
Conversely,  a  parameter  to  which  the  model  is  highly  sensitive  may  contribute 
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substantially  to  the  output  uncertainty  even  if  the  uncertainty  associated 
with  the  parameter  is  relatively  minor. 

B .  METHODS 


The  technique  that  was  used  to  measure  the  contribution  of  individual 
parameters  involved  estimation  of  the  variance  of  the  output  (dose  surrogate 

or  risk  estimates).  The  variance  for  a  sample  of  N  values  (x^,  X2 .  x„) 

is  defined  as 


(1) 

V  -  l(x,  -  p)/(N-l) 

C2) 

H  -  Ixi/N 

where  the  suras  in  each  equation  are  over  all  values  of  i  (-1,  2 . N)  . 

Let  denote  the  variance  associated  with  the  output  variable  values 
when  all  parameters  were  allowed  to  vary  according  to  their  uncertainty 
distributions.  Similarly,  let  denote  the  variance  associated  with  the 
output  variable  values  when  parameter  i  (or  parameter  set  i)  was  allowed  to 
vary,  but  all  other  parameters  were  fixed.  Then 

(3)  Ri  -  Vi/Vj 


is  a  representation  of  the  contribution  of  parameter  i  (or  parameter  set  i) 
the  overall  uncertainty  associated  with  the  output  variable  values.  When  Rj 
is  small,  parameter  i  has  a  small  Impact  on  uncertainty  In  the  output 


to 
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variable.  Larger  values  of  indicate  greater  impact  on  the  variability  of 
the  output  variable. 

The  representation  character izad  by  the  valu  s  combined  two  factors. 
One  was  the  amount  of  variation  allow^'d  for  a  parameter,  as  represented  by  its 
uncertainty  distribution.  Clearly,  if  a  parameter  was  assumed  to  be  known 
with  certainty  (i.e.,  the  uncertainty  distribution  allowed  for  no 
variability) ,  then  R^  for  that  parameter  would  be  zero  (since  would  be 
zero).  The  more  uncertain  a  parameter  value  is  (i.e.,  the  greater  the 
variation  allowed  by  its  associated  uncertainty  distribution)  the  greater  the 
chance  that  the  R,  for  that  parameter  will  be  larger. 

However,  the  other  factor  that  affected  the  values  of  R^  influenced  the 
degree  to  which  increasing  parameter  uncertainty  is  manifested  in  increased 
output  variable  uncertainty.  The  second  factor  was  the  sensitivity  of  the 
model  to  the  parameter  under  consideration.  Thus,  a  parameter  to  which  the 
PBPK  model  is  relatively  insensitive  will  probably  have  a  small  R^  value,  even 
if  the  uncertainty  associated  with  that  parameter  is  very  large.  Sensitivity 
for  the  perchloroethylene  PBPK  model  was  addressed  in  detail  in  Part  1. 

C .  RESULTS 

Tables  VI-2-1  and  VI-2-2  display  results  that  extend  the  analysis  of  the 
perchloroethylene  example  presented  in  Appendix  Vl-l-A.  Those  tables  show  the 
Rj  values  for  each  parameter  or  set  of  parameters  when  the  output  variable  was 
either  a  dose  surrogate  estimate  or  a  risk  estimate.  The  dose  surrogates 
examined  include  those  discussed  in  Part  1  of  this  volume  and  Appendix  VI-l-A. 
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In  one  instance  (the  contribution  of  uncertainty  In  pla  to  uncertainty 
in  the  estimation  of  human  AUCL;  Table  VI-2'1),  the  variability  in  the  output 
found  when  a  single  parameter  was  varied  exceeded  that  observed  when  all 
parameters  varied.  This  result  may  have  arisen  as  follows.  When  all 
parameters  were  allowed  to  vary,  the  random  matching  of  sampled  parameter 
values  may  have  moderated  some  of  the  extreme  AUCL  values  produced  when  pla 
alone  varied.  That  would  tend  to  decrease  the  variability:  l.e.,  AUCL 
variance  associated  with  uncertainty  in  the  pla  parameter  would  be  smaller 
than  that  associated  with  uncertainty  in  all  the  parameters,  as  observed. 

TV'i  Q  will  examined  further. 

The  contributions  of  all  of  the  individual  parameters  or  sets  of 
parameters  to  risk  uncertainty  were  small  when  virtual  concentration  of 
metabolite  was  used  as  the  dose  surrogate  (Table  VI-2-2).  This  can  be  seen  by 
comparing  values  in  the  three  columns  of  Table  VI -2-2.  This  result  was 
accompanied  by  a  shift  in  the  median  risk  estimates;  i.e.,  the  median  risk 
estimate  when  all  parameters  were  allowed  to  vary  was  higher  than  the  risk 
estimate  medians  observed  when  individual  parameters  were  allowed  to  vary. 

This  result  could  have  arisen  in  the  following  way.  The  original  fit  of 
the  multistage  model  (using  the  observed  response  rates  and  the  preferred  PBPK 
parameter  values)  estimated  that  the  linear  term  of  Che  model  was  zero. 

Changes  in  parameter  values  affect  risk  estimates  by  altering  multistage  model 
parameter  estimates,  but  such  changes  could  never  decrease  the  linear  term. 
Consequently,  risk  estimates  would  tend  to  be  increased,  especially  at  low 
doses  where  the  linear  term  dominates.  Thus,  when  all  parameters  were  allowed 
to  vary,  the  linear  term  tended  t''  be  increased  more  often  than  not  because  of 
the  overall  variability.  However,  when  individual  parameters  were  varied,  the 
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variability  of  the  dose  surrogate  estimates  was  not  so  great,  inducing  a 
smaller  change  in  risk  estimates  (l.e.,  tending  not  to  increase  the  linear 
term  more  often  than  not) ,  This  would  cause  a  shift  in  the  median  risk 
estimates  (smaller  median  values,  more  similar  to  that  obtained  with  the 
original  fit  of  the  multistage  model)  when  single  parameters  varied,  and  may 
also  have  resulted  in  less  variability  associated  with  the  risk  estimates. 

A  similar  phenom.enon  was  not  observed  and  would  not  be  expected  for  risk 
estimates  based  on  the  other  two  dose  surrogates.  This  is  because,  in  those 
cases,  the  linear  term  was  positive  when  the  multistage  model  was  fit  to  the 
observed  responses  using  the  preferred  PBPK  model  parameters.  Thus,  changes 
in  the  parameters  (either  one  at  a  time  or  all  together)  could  induce  both 
increases  or  decreases  in  the  linear  term,  tending  to  keep  the  median  close  to 
that  observed  with  the  original  fit. 

Despite  the  fact  that  the  above  arguments  can  explain  the  pattern  of 
risk  estimates  observed  for  risks  based  on  the  metabolite  vir.ual 
concentration,  further  investigation  of  this  issue  will  be  carried  out. 

D.  DISCUSSION 

The  contributions  of  the  individual  parameters  to  the  output 
uncertainties  matched  fairly  closely  the  pattern  observed  when  the  model  was 
examined  with  respect  to  sensitivity  (Part  1  of  this  volume).  Contribution  to 
dose  surrogate  uncertainty  corresponded  extremely  well  with  contribution  to 
risk  estimate  uncertainty. 

For  the  AUCL  dose  surrogate,  the  model  was  particularly  sensitive  to  pla 
and  pb .  Those  two  parameters  contributed  meet  highly  to  the  AUCL  dose 
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surrogate  uncertainty  (Table  VI -2-1).  This  was  true  for  the  human  pb  and  pla 
values  as  well  as  for  the  mouse  pb  and  pla  values. 

For  the  dose  surrogate  based  on  virtual  concentration  of  metabolite,  the 
response  rates  contributed  most  highly  to  the  risk  uncertainty.  The  metabolic 
constants  were  the  most  significant  contributors  to  risk  uncertainty  of  all 
the  PBPK  model  parameters.  Because  of  the  relatively  low  atmospheric 
concentrations  to  which  the  mice  were  exposed  (100  and  200  ppm),  the  mouse  km 
parameter  contributed  more  to  risk  uncertainty  than  did  the  mouse  vmaxc 
parameter.  Model  predictions  of  virtual  concentration  of  metabolite  were 
found  to  be  very  sensitive  to  the  value  of  the  liver  volume  parameter  (Part  1 
of  this  volume) .  That  sensitivity  is  reflected  in  the  relatively  large 
contribution  of  uncertainty  in  all  volumes  (varied  together  because  of  their 
interrelationships)  to  dose  surrogate  and  risk  uncertainty,  when  based  on 
metabolite  virtual  concentration. 

The  dose  surrogate  AUCA  was  most  sensitive  to  the  blood/air  partition 
coefficient.  This  parameter  also  contributed  very  significantly  to 
uncertainty  in  the  risk  estimates.  The  contributions  of  uncertainties  with 
respect  to  response  rates,  mouse  vmaxc,  and  mouse  km  were  also  substantial  for 
the  risks  based  on  the  AUCA  dose  surrogate.  The  contributions  of  all  other 
parameters  were  considerably  less. 

Some  parameters  were  consistently  unimportant  in  determining  output 
uncertainties.  Those  parameters  included  the  partition  coefficients  for  the 
fat,  slowly  perfused,  and  richly  perfused  compartments.  Of  those  partition 
coefficients,  the  one  for  the  slowly  perfused  compartment  contributed  to 
output  uncertainty  more  than  the  other  two.  The  parameter  describing  the 
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relationship  between  the  waking  and  sleeping  cardiac  output  rates  In  humans 
(dqc)  was  also  consistently  unimportant  for  determining  output  uncertainty. 
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Table  Vl-2'1 


Contributions  of  Individual  Parameters 
to  Uncertainty  in  Dose  Surrogate  Estimates* 


_ Stirrogate^ _ 

Virtual  Concentration 


Parameter 

AUCL 

of  Metabolite 

AUCA 

Vmaxc 

u 

6E-2 

(^) 

19 

(2) 

4. 

7E-1 

(3) 

Km 

f. 

.OE-2 

(3) 

22 

(1) 

7  . 

OE-1 

(2) 

pb 

3 

.  8 

(2) 

4  . 

1 

(3) 

97 

(1) 

p  1  a 

1 

X. 

22 

(1) 

5. 

OE-2 

!  (7) 

7. 

OE-2 

(6) 

pra 

.  S  E  -  3 

(") 

5 

1E-: 

>  (6) 

6  . 

8E-2 

(8) 

psa 

1 

.  lE-2 

(5) 

4  . 

8E-  ] 

.  (5) 

3. 

4E-1 

(9) 

p  f  a 

1 

.  7E-3 

(8) 

4  . 

OE-1 

!  (9) 

6  , 

9E-2 

(2) 

Volumes'^ 

1 

.  5E-3 

(9) 

1  . 

7 

(9) 

6  . 

lE-2 

(9) 

doc"' 

0 

.  t  L-  3 

(6) 

4  . 

6E-i 

’  (8) 

8  . 

lE-2 

(5) 

’  Presented  are  values  of  Ri*100  (see  text  for  definition)  and,  in 

parentheses,  a  rank  for  the  parameter.  A  rank  of  1  indicates  largest  Ri . 
Tliese  are  for  the  human  PF.RC  model  with  a  50  ppm,  8  hour/day,  5  day/week 
e  '  sure  .sc  ena  r  i  o 

“  The  dose  surrogates  are  average  daily  area  under  the  liver  concentration 
curve  (AUCL) ,  average  daily  amount  metabolized  per  liver  volume  (virtual 
concentration  of  metabolite),  and  average  daily  area  under  the  arterial 
blood  concentration  curv’e  (AUCA)  . 

'  The  volumes  are  correlated  (they  vary  together  according  to  a  Dirichlet 
distribution)  and  so  tl,f-ir  cunt  r ;  hut  ion  to  cutput  uncertainty  constitutes 
one  entry. 

The  parameter  defining  the  relationship  between  resting  and  active  cardiac 
output  rates.  See  Part  1  of  this  volume. 
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TaV^le  VI -2-2 


Contributions  of  Individual  Parameters 
to  Uncertainty  In  Risk  Estimates* 


_ Surrogate*’ _ 

Virtual  Concentration 

Pa raireter  _ AUCL _ of  Metabolite _ AUCA 


Human 

Vmaxc 

2  . 

OE- 

1 

(10) 

1 , 

3E 

3 

(3) 

1. 

3E- 

•1 

(7) 

Hvumn 

Kjn 

2 

6E  - 

1 

(9) 

6. 

OE 

4 

(5) 

1. 

9E- 

•1 

(6) 

Human 

pb 

13 

(2) 

4 

6E 

5 

(6) 

32 

(1) 

Human 

pi  a 

88 

(1) 

3 

3E 

7 

(12) 

2  . 

OE- 

■2 

(15) 

fiunian 

p  r  a 

7  . 

5E- 

3 

(16) 

3, 

2E 

7 

(13) 

1. 

9E- 

■2 

(17) 

Human 

psa 

4  . 

6E- 

2 

(14) 

3  . 

8E 

6 

(10) 

9  . 

7E- 

•2 

(13) 

Human 

pf  a 

7  . 

4E- 

■3 

(17) 

2 

5E 

7 

(16) 

1 

9E 

2 

(16) 

H’.mian 

volume' 

6  . 

lE- 

3 

(18) 

1 , 

5E 

5 

(7) 

1. 

7E- 

•2 

(18) 

Human 

dqc'* 

1  . 

OE- 

■2 

(15) 

2, 

9E 

7 

(14) 

2. 

3E- 

•2 

(14) 

.^ouse 

Vmaxc 

8  . 

8 

(5) 

7 

7E 

4 

(4) 

9  . 

1 

(4) 

.Mouse 

Km 

X  . 

0 

(7) 

1 

7E 

2 

(2) 

2. 

3 

(5) 

Mouse 

pb 

13 

(3) 

5 

,2E 

6 

(8) 

26 

(2) 

Mouse 

pla 

13 

(4) 

1 

,4E 

8 

(18) 

1 , 

OE- 

■1 

(10) 

Mouse 

p  r  a 

7 . 

6E- 

•2 

(13) 

1 

.9E 

8 

(17) 

1. 

OE- 

•1 

(10) 

Mouse 

psa 

9. 

9E- 

•1 

(8) 

2 

.8E 

7 

(15) 

1 . 

2E- 

•1 

(8) 

Mouse 

pfa 

1. 

OE- 

•1 

(11) 

6 

.9E 

7 

(11) 

1. 

lE- 

■1 

(9) 

Mouse 

volumes' 

8  . 

OE- 

■2 

(12) 

3 

.9E- 

■6 

(9) 

1 . 

OE- 

■1 

(10) 

Response  rates® 

3. 

0 

(6) 

5 

.0 

(1) 

11 

(3) 

*  Presented  are  values  of  R^*100  (see  text  for  definition)  and,  in 
parentheses,  a  rank  for  the  parameter.  A  rank  of  1  indicates  the  largest 
Rj .  Risks  were  derived  for  a  50  ppm,  8  hour/day,  5  day/week  human  exposure 
scenario  using  the  female  mouse  results  discussed  in  Appendix  VI-l-A. 

The  dose  surrogates  are  average  daily  area  vr.^er  the  liver  concentration 
curve  (AUCL) ,  average  daily  amount  metabolized  per  liver  volume  (virtual 
concentration  of  metabolite),  and  average  daily  area  under  the  arterial 
blood  concentration  curve  (AUCA) . 

The  volumes  are  correlated  (they  vary  together  according  to  a  Dirichlet 
distribution)  and  so  their  contribution  to  output  uncertainty  constitutes 
one  entry . 

The  parameter  defining  the  relationship  between  resting  and  active  cardiac 
output  rates.  See  Part  1  of  this  volume. 

*  The  bioassay  response  rat^s  were  allowed  to  vary  according  to  a  binomial 
distribution  with  probability  of  response  defined  by  the  multistage  model 
fit  to  observed  response  rates  with  doses  defined  by  preferred  values  of 
PBPK  model  parameters. 
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